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Abstract 

Discussed is kinematics and dynamics of bodies with affine degrees of 
freedom, i.e., homogeneously deformable "gyroscopes". The special stress 
is laid on the status and physical justification of affine dynamical invari- 
ance. On the basis of classical Hamiltonian formalism the Schroedinger 
quantization procedure is performed. Some methods of the partial separa- 
tion of variables, analytical treatment and search of rigorous solutions are 
developed. The possiblity of applications in theory of structured media, 
nanophysics, and molecular physics is discussed. 

Keywords: affine degrees of freedom, molecular dynamics, nanophysics, quan- 
tized media, structured media. 

Introduction 

The mechanics of affine bodies was a subject of many papers [SI EH E2 Efl EE 

ng nzi hhi noi mj usi na us eoi esi ezi ehi esi sh ssi sn S3 sa 

US [Ml ISH [331 [511 1531 [351 [371 [BOl EH 153 EH US JES [H71 Ha [BSl [701 EH 

[721 Ea EH US [7B1 E3 Ea HOI- It has been a field of intensive studies in our 
group at the Institute of Fundamental Technological Research in Warsaw. Up 
to our knowledge, for the first time the idea of objects with affine degrees of 
freedom in mechanics appeared in papers of Eringen [26l [27j [28l [29] devoted to 
structured continua, to be more precise in his theory of micromorphic media. 
Micromorphic continuum is an affine extension of the micropolar Cosserat con- 
tinuum. Roughly speaking, the Cosserat medium is a deformable continuum 
of infinitesimal gyroscopes. Similarly, the micromorphic body is a deformable 
continuum of infinitesimal homogeneously deformable gyroscopes. Affine model 
of collective degrees of freedom was also used in the theory of collective phenom- 
ena in atomic nuclei Ej . The idea of affine body is interesting in itself from the 
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point of view of analytical mechanics and theory of dynamical systems. It is an 
instructive example of systems with degrees of freedom ruled by Lie groups. In 
mechanics of non-constrained continua the configuration space may be identified 
with the group of all diffeomorphisms of the physical space (volume-preserving 
diffcomorphisms in mechanics of ideal incompressible fluids). It is rather diffi- 
cult to be rigorous with such infinite-dimensional groups. Affine model is placed 
between rigid-body mechanics and the general theory of deformable continua, 
i.e., it involves deformations but at the same time one deals there with a finite 
number of degrees of freedom. The Lie-group background of the geometry of the 
configuration space offers the possibility of the effective use of powerful analytic 
techniques. One can realize certain finite-dimensional generalizations when the 
configuration space geometry is ruled, e.g., by the projective or conformal group. 
Also other finite-dimensional discretized approaches are useful but of course the 
models based on geometric transformation groups are particularly interesting 
and efficient. 

The range of applications of affine model of collective and internal degrees of 
freedom is very wide and has to do with various scales of physical phenomena: 

• macroscopic elastic problems when the length of excited waves is compa- 
rable with the linear size of the body. 

• purely computational and engineering problems connected with the finite 
elements methods. A mixture of analytic and numerical procedures. 

• structured bodies, e.g., micromorphic continua and molecular crystals. 

• vibrations of astrophysical objects (stars, concentrations of the cosmic 
dust), theory of the shape of the Earth. 

• molecular vibrations. 

• nuclear dynamics. 

Obviously, the last two subjects must be based on the quantized version of 
the theory. Quantum description is also necessary in various problems concern- 
ing the nanoscale phenomena, fullerens, etc. It is a new fascinating subject 
where one deals with the very intriguing convolution of the classical and quan- 
tum levels, perhaps also with some yet non-solved paradoxes from the realm of 
quantum-mechanical foundations like decoherence, etc. 

Quantization as a purely mathematical procedure is connected with certain 
ambiguities which may be solved only a posteriori, on the bases of experimental 
data. There arc some well-known problems with the ordering of operators. 
In models with a firm group-theoretic background there are some canonical 
procedures, usually confirmed by experiments. Because of this, the extensive 
geometric introduction presented below, almost a treatise as a matter of fact, 
is a constitutive element of the theory, motivated by deeper reasons than the 
purely mathematical curiosity or artificial sophistication. 

Affine models of degrees of freedom of the structured elements is very natural. 
When one deals with fullerens, macromolecules, microdefects, affine modes of 
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motion are certainly the most relevant ones. There are molecules, e.g., P4, which 
have no other degrees of freedom; there are also such ones for which non-affinc 
behaviour is a merely small correction. As mentioned, affine modes have also 
to do with finite elements, when the body is described as an aggregate of small 
affine objects. 

Perhaps the quantization of such an approach might be a procedure alter- 
native to the phonon description based on the quantized plane elastic waves. 

And finally, one of the most important things. The group-theoretical de- 
scription of internal and collective modes is really effective when the dynamics 
is invariant (or in some sense almost invariant) with respect to the group under- 
lying kinematics of the problem. And this is not the case in all models of affine 
bodies met in literature. Kinematics is there affine but the group of dynamical 
symmetries is broken to the isometry group. Because of this, there is no full 
use and the full profit of Lie-group techniques. Unlike this, we formulate here 
affinely-invariant dynamics, where elastic interactions may be encoded in appro- 
priate kinetic energy models without (or "almost" without) any use of potential 
energy terms. This procedure is similar to that following from the Maupertuis 
variational principle. There are indications that just such models may be useful 
in condensed matter theory, where the structural elements are more sensitive 
to the geometry of a surrounding piece of the body, e.g., to the Cauchy defor- 
mation tensor than to the "true" metric tensor of the physical space. This is 
something similar to the effective mass tensors of electrons in crystals. 

There is an interesting link between our models and theories of integrable 
lattices like Calogero-Moser, Sutherland, and others [40l [41] . On the quan- 
tum level the deformation invariants behave like indistinguishable, exotically 
parastatistical one-dimensional " particles" . 

Obviously, the real world, the arena of mechanical phenomena, is three- 
dimensional. However, certain important invariance and other problems are 
explained in a more lucid way when described with non-physical generality, i.e., 
in n dimensions. By the way, two-dimensional problems are also interesting 
not only in "Flatland" [T] but also in some realistic physical problems. At 
the same time, they are computationally simple due to some exceptional, so to 
speak pathological, feature of GL(2,R) among all GL(n,R) (SO(2,R) is Abelian, 
whereas SO(n,R) for n > 2 are semisimple). 

1 Classical preliminaries 

Let us briefly describe various models of the configuration space of affinely rigid 
body, ft depends on the particular problem under consideration which of them 
is more convenient. The possibility and usefulness of many choices of geometric 
structures underlying physically the same degrees of freedom was pointed out 
by Capriz [9j [lOl QTJ [12l [33] [34j [35] . Various descriptions differ in assuming 
some auxiliary geometric objects. 

We begin with some elementary concepts of affine geometry, just to fix the 
language and notation. Affine space is given by a triple (A, E, — >), where X is 
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a point set, just the "space itself, E is a linear space of translations in X, and 
the arrow — > denotes a mapping from the Cartesian product X x X onto E; 
the vector assigned to (p, q) G X x X is denoted by pH}. The arrow operation 
satisfies some axioms, namely, 

(i) p$ + qf + rf> = for any p, g,r£ X, 

(ii) for any p G X and v G E there exists exactly one q £ X such that pq~ = v; 
we write q = t v (p) . 

For any v G E, t v : X — > X is a one-to-one mapping of X onto X, the translation 
by v. And obviously, 

t v o t u t u o t v t v -^- U: t Q — idx , t u t—u- 

In this way, E considered as an additive-rule Abelian group acts freely and 
transitively on X. 

Any linear space E may be considered as an affine space (E,E,—), i.e., 
ui — v — u. 

The axiom (i) implies that pf> — 0, pq" — —qp for any p,q G X. 

Let f2 be an arbitrary set, in general, structure-less one. The set of all 
mappings from f2 in X, denoted by X^, is simply the 17-indexed Cartesian 
product of X. For any mapping / : O — > X, the image f(ui) e X is interpreted 
as an w-th component of /. When is a finite A^-element set, e.g., il = 
{1, 2, . . . , AT}, this is just the familiar finite Cartesian product X N . 

The set X n is in a natural way an affine space. Its translation space is 
identical with E n , the set of all mappings from SI into X. I f F l G are mappings 
from Q into X, then the corresponding translation vector FG G is simply 
given by 

(FG) ( W ) := F(c)G(cj, (1) 

for any u; € O. One can easily show that all axioms of affine geometry are 
satisfied then. 

Affine mappings, by definition, preserve all affine relationships between fig- 
ures and points. So, if (A, U, — (M, V", — >) arc affine spaces, then we say that 
$ : A — > M is affine if there exists such a linear mapping L[<f>] : U — > V denoted 
also by that for any p,g£ JV the following holds: 

*(p)*(gj = i[*]p3- 

The mapping L[<f>] : [/ — > V is referred to as a linear part of $. The set of all 
affine mappings from N to M will be denoted by Af(A, M); similarly, L(£7, V) 
denotes the set of linear mappings. If $i G Af(P, M) and $2 G Af(A, P), then 
$1 o $ 2 g Af(A, M) and o $ 2 ] = i[$i]L[$ 2 ]- 

Dimension of the translation space E is referred to as the dimension of X 
itself. Any fixed point p G M establishes the bijection of M onto V given 
by Mb q h € V. Such V- valued charts establish in M the structure of 
analytical differential manifold just of dimension dim V. The manifold of affine 
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injections from N into Af will be denoted by AfI(iV, Af), and the corresponding 
set of linear injections from U into V by LI(C/, V). They are open submanifolds 
of Af(N,M), L(U,V), respectively. Obviously, they are non-empty if and only 
if dim Af > dim N. If dim Af — dim N, they become respectively the manifolds 
of affine and linear isomorphisms. 

If N = M and U — V, i.e., when we work within some fixed affine space 
(Af, V, — »), then some simplified notation is used, namely, 

L(V, V), Af (Af, M), U.(V,V), AR(M,M) 

are denoted respectively by 

L(V), Af (M) , GL(V), GAf(Af). 

Obviously, the last two sets are groups, respectively, the general linear and affine 
groups in V, M. Translations are affine isomorphisms; their set T[V] = {t v : 
v G V} is a normal subgroup of GAf(M). This subgroup is the kernel of the 
group epimorphism: 

GAf (Af) 3^ L[cp] e GL(V). 

The quotient group GAf(M)/T(V) is isomorphic with GL(V) but in a non- 
canonical way; any choice of centre o G Af gives rise to some isomorphism. 

The set of affine mappings from (N, U, — ►) to (Af, V,—*), i.e., M(N,M), is 
an affine subspace of M N in the sense of (JXJ) ; the translation space is identified 
with V N . 

If in the space N some origin point 6 TV is chosen then the manifold 
Ai(N,M) may be simplified to the Cartesian product Af x L(U,V). Namely, 
with any $ 6 Af(iV, M) we associate a pair (x, <y?) G Af x L(C7, V) in such a way 
that x = <I>(0) and $(C) ( I > (a) = 99 • Oa. When we restrict ourselves to the open 
submanifold of affine isomorphisms AfI(7V, Af) C Af(iV, Af ), then ip in the above 
expression runs over the open submanifold Ll(U, V) C L(t/, V). 

And finally, let us fix some linear frame, i.e., an ordered basis in U, E = 
(E±, . . . , Ea, ■ ■ ■ , E n ), n = dim U — dim V. When it is kept fixed, any linear 
mapping ip G h(U, V) may be identified with the system e = (ei, . . . , eji, . . . , e„), 
where e A = V>Ea, A = T~n. When $ G AfI(AT, Af), i.e., ip G L1(U, V), then e is 
a linear frame in V . In this way LI(U, V) is identified with F(V), the manifold 
of linear frames in V. And AfI(A, Af ) is identified with Af x F(V), the manifold 
of affine frames in Af (the pairs consisting of points in Af and ordered bases in 
V). 

Fixing an affine frame (O, E) in N we turn it into the arithmetic space 1". 
Linear isomorphisms of U onto V become then linear frames in V; their inverse 
isomorphisms are identified with the dual co-frames: e — (e 1 , . . . , e A , . . . , e n ), 
(e A ,e B ) = 5 A B 

As frames and dual co-frames mutually determine each other, AfI(A r , Af ) 
may be as well identified with Af x F(V*) = Mx F(V)*; here F(V*) is the 
manifold of frames in the dual space V* denoted also as F(V)*. 
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If in addition some affine frame (o, £) — (o; £\, . . . , £a 7 ■ ■ ■ , £ n ) in M is fixed, 
then also V becomes identified with R™. The manifold Ll(U, V) is then analyt- 
ically identified with the general linear group GL(n, M), and AfI(iV, M) may be 
identified with the semi-direct product GAf(n,K) ~ GL(n,R) x s R n . 

The manifold AfI(iV, M) is a homogeneous space of affine groups GAf(M), 
GAf(AT) acting respectively on the left and on the right: 

A £ GAf(M) : AfI(A, M) 3$hAo $, (2) 
B £ GAf (AT) : AfI(A, M)9$h$oB. (3) 

Similarly, linear groups GL(V), GL(/7) act transitively on Ll(U, V): 

a £ GL(V) : IA{U, V) 3 cp h-> aip, (4) 
G GL(J7) : LI(C/, V) 3 ip i-> (5) 

Let us observe that although GL(V), GL({7) are logically distinct disjoint 
sets, the corresponding transformation groups intersect nontrivially. Namely, 
dilatations belong to both of them, the left and right actions of a = X Idy, 
(3 = A Id[/ result in multiplying (p by A, i.e., (p > Xip. 

When some origin O e JV is fixed and Afl( N, M) is identified with M x 
LI(U,V), the left-acting transformation groups GAf(M), GL(F) may be repre- 
sented as follows: 

A e GAf (M) : M x LI(U, V) 3 (x, (p) ^ (A{x), L(A)<p), (6) 
a 6 GL(V) : M x Ll(U, V) 3 (x, <p) ^ (x, aip), (7) 

The origin O enables one to identify GAf(iV) with the semi-direct product 
GL(U) x s U. Namely, B S GAf(JV) is represented by the pair (L(B), OB(0)). 
And conversely the pair (fl, u) G GL(U) x s U gives rise to the mapping B G 
GAf(iV) such that for any a £ N: 

OB (a) =/3-Oa + u. 

The right action §5§ of B on AfI(A, M) is represented in M x LI(U, V) as follows: 

{j3,u)£GL{U)x s U : (a:, ^) » (t vu )(x), <pp). (8) 

If we put u = 0, then the group GL([/) itself acts only on the second component: 

p£GL(U) : (x,<p)^(x,<p(3). (9) 

The standard language of continuum mechanics is based on the use of two 
affine spaces: the physical and material ones. We denote them respectively by 
(M, V, — >) and (N, U, — >). If we deal with the infinite continuum medium filling 
up the whole physical space, configurations are described by diffeomorphisms 
of N onto M. The smoothness class of these diffeomorphisms depends on pe- 
culiarities of the considered problem. The manifold N is interpreted as the set 
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of material points. In configuration given by $ : -/V — > M, the material point 
a e N occupies the spatial position $(a) € M. Diffcomorphism groups Diff(M) 
and Diff(iV) give rise to transformation groups acting on the configuration space 
Diff(JV,M), i.e., 

AeBiS(M) : Diff (N, M) 3 $h4o$, (10) 
BeVitt(N) : Diff (TV, M) 3 $ $ o B. (II) 

They are referred to respectively as spatial and material transformations. Ob- 
viously, spatial and material transformations mutually commute. In continuum 
mechanics they have to do with symmetries of space and material itself. 

Obviously, when one deals with realistic bounded bodies, this description 
should be modified, e.g., manifolds with boundary become a better model of 
the material space. Another possibility is to use a smooth smeared-out model 
of the boundary, i.e., to describe the bounded body as a non-bounded, one 
however, with the mass density quickly vanishing outside the real object. 

Deeper modifications are necessary when describing continua with degen- 
erate dimension like membranes, strings, infinitesimally thin shells, rods, etc. 
And obviously, for discrete systems the description based on the affine space N 
as a material body is not applicable in the literal sense, unless some tricks like 
smeared out density functions and so on are used. It would be a good thing, 
especially when dealing with affine systems in microscopic applications (molec- 
ular, microstructural, etc.) to start from some general formulation applicable 
both to discrete and continuous systems of various kinds. There is also some 
more subtle point to that, namely, the material space is primarily the abstract 
set of material points or their labels, so-to-say "identification cards". A priori 
this set is structure- less; it is a kind of "powder" of material points. Let us 
denote it by 0. Configurations are mappings from Q, to M, i.e., elements of 
M n (the usual finite Cartesian product M N when one deals with an iV-particle 
system). More precisely, M n is the set of singular configurations, i.e., ones ad- 
mitting coincidences of different material points at the same spatial point. To 
avoid such a "catastrophe" one must decide that the "true" configuration space 
is the set of injections from into M, i.e., Inj(f2, M). As far as Q, is structure- 
less the only well-defined set of material transformations is Bij(il), the set of 
bijections of f2 onto f2. They act on configurations according to the following 
rule: 

B e Bij(fi) : Inj(fi, M) 3 $ ^ $ o B. (12) 

These "permutations" of material points are just the only admissible material 
transformations on this yet amorphous stage. Transformations of M onto itself, 
in particular diffeomorphisms of M onto itself of an appropriate class, act on 
Inj(f2, M) according to the following rule: 

AeBij(M) : Inj(f2, M) 9$h4o$. (13) 

In general, this action is not transitive and splits into orbits, i.e., transitively 
classes. Any fixed class carries over geometric structures from M to O. For 
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example, if f2 has the continuum cardinal number and only the bijections of f2 
onto M are admitted as configurations, then any fixed orbit of the left-acting 
diffeomorphism group DifF (M) induces in £1 some structure of C r -class differ- 
entiable manifold. The powder of material points becomes the continuous body 
and its configuration space is identified with DifF (fl,M), i.e., the set of C r - 
class diffeomorphisms of il onto M. If some orbit of the left-acting affine group 
GAf(Af) is fixed, f2 becomes endowed with the induced structure of affine space. 
And the one can sensibly tell about affine mappings from f2 onto M and about 
affine relationships between material points in Q,. 

Different orbits induce structures in Q, which are literally different, although 
usually isomorphic. 

In general, when dealing with constrained systems of material points, it is not 
the total group Bij(M) or the diffeomorphism group DifF (M), but some rather 
peculiar proper subgroup G C Bij(M) that rules geometry of degrees of freedom 
and perhaps also the dynamics. Configuration spaces are constructed as orbits 
of such groups. Let us assume that some orbit Q, i.e., some particular model of 
degrees of freedom is fixed. By the very definition, G acts transitively on Q, i.e., 
Q is a homogeneous space of the action (fl3]) . The point is how to define some 
right-hand-side action analogous to (|II[) or In general, the transformation 
group Bij(f2) acting through (Tl^|) is too poor. For example, when one deals 
with a finite system of material points, Bij(il) is the permutation group of f2 
and nothing like continuous groups of material transformations, e.g., (|f I[) . ^ 
can be constructed on the basis of Bij(f2). One feels intuitively that there is 
something non-satisfactory here. And indeed, it is possible to define some rich, 
in general continuous, group of material transformations acting on the right on 

Q- 

The construction is more lucid when one forgets for a moment about details 
and considers an abstract homogeneous space with the underlying point set Q 
and the group G acting transitively on Q on the left. For simplicity, we denote 
the action of G by 

g E G : Q 3 q i— > gq E Q. 

This graphical convention is well-suited to the left-hand-side nature of this ac- 
tion: 

{gi92)q = 91(920)- 

The action is assumed to be effective, i.e., the group identity e E G (its neutral 
element) is the only element of G satisfying the condition eq = q for any q E 
Q. Roughly speaking, G is a "proper" transformation group, not something 
homomorphically (with a non-trivial kernel) mapped into the group Bij(Q). 
This is exactly the case in problems described previously. Let qo E Q denote 
some arbitrarily fixed points, and H(qo) C G denote its isotropy group, i.e., the 
set of elements which do not move q: 

H(q ) := {g EG : gq Q = q } . 
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It is well-known that Q may be identified in a one-to-one way with the set of 
left cosets, i.e., with the quotient space G/H(q ). The action of G on G/H(q ) is 
represented by left translations, namely, the coset xH(q n ) := {xh : h <G H(q )} 
is transformed by g € G into gxH(qo); obviously, the result does not depend 
on the particular choice of x within its coset, i.e., on the replacement x i— > xh, 
heH. 

The question is now whether there exist some right translations of represen- 
tants, x i ► xg, admitting an interpretation in terms of transformations acting 
on G/H(q ). It is easy to see that the answer is affirmative. 

Namely, let N (H(qo)) C G denote the maximal subgroup of G for which 
H(q a ) C N (H(q )) is a normal subgroup. It is easy to see that for every 
n e N (H(q j) the corresponding right regular translation G 3 x xn e 
G is projectable to the manifold of left cosets G/H(q ). Indeed, cosets are 
transformed onto cosets: 

(xH(qo)) n — xH(qo)n = x (nH(qo)n^ 1 ^ n = xnH(qo) = (xn)H(qo). 

In this way transformations may be performed on representants, G 3 x xn e 
G. Obviously, the choice of representants does not matter because 

(xh)nH (qo) = (xhn)H(qo) ~ xnn^ 1 hnH (qo) , 

and for any h e H(qo), n^ 1 hn S H(qo), thus n~ 1 hnH(qQ) = H(qo), and finally 
((xh)n)H(qo) = (xn)H(q ). The non-effectiveness kernel of the right action of 
N (H(q )) on the coset manifold coincides with the group H(q ) itself, thus the 
true group of right-acting transformations is given by 

N (H(q )) /H(q ) = H(q ) \ N (H(q )) . 

The above construction of right-acting transformations pre-assumes some choice 
of the reference point go € Q- The question arises as to what extent does the 
presented prescription depend on q . It turns out that the constructed transfor- 
mation group itself is well-defined and the particular choice of qo influence only 
the "parameterization", so-to-say identification labels of the group elements. 
Let gi,<72 G Q be two arbitrarily chosen reference point. The subset of G con- 
sisting of elements k transforming qi into 92, i-e., kq\ = q2, will be denoted by 
H{qi,q2)- Obviously, H(qi,q 2 ) is simultaneously the left and right coset of the 
subgroups H(qi), H (92), respectively; if k is an element of H (51,52), then so 
is hikh\ for any h\ G H(qi), h 2 € H(q 2 ). (Incidently, let us notice: perhaps 
it would be convenient to write H(qi,qi), H(q 2 ,q2) instead of H(qi), H{q 2 ).) 
Obviously, for any k e H(qi,q 2 ) we have 

H(q 2 ) = kH(q 2 )k-\ N(q 2 ) = kN^)^ 1 . 

Any choice of k e H(q\,q 2 ) fixes some isomorphisms of H{q{), N(q\) respec- 
tively onto H(q 2 ), N(q 2 ). 

Let us take some g\ e G and the point q e Q produced by it from qi e Q, 

i.e., 

q = 9i1i- 
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When q is fixed, g\ is denned up to the gauging g\ i— > g\h, h 6 H(qi). And now 
we transform g by the right action of n\ € iV (H(qi)): 

qh^q' = gimqt. 

The result does not depend on the gaugings <?i i— ► tfa/i, ri\ i— > Xi^iXi, where 
^> Xi) X2 € H(qi). Let us now express this action in terms of the reference point 
qi = kqi 

q' = giniqi = g\n 1 k~ x q 2 = (gik^ 1 ) (knik^ 1 ) q 2 . 

Now q is produced from the reference point q 2 by g 2 — gik^ 1 ; q — g 2 q 2 . And 
its representing group elements g 2 is affected on the right by n 2 — kn\k~ x € 
N (H(q 2 )) the fc-conjugation of n\ € N(H(qi)). In this way different right- 
hand-side actions in G, i.e., 

Gsih xm, Gain xn 2l 

describe the same transformation in Q. They are different "labels" of this 
transformation corresponding to various choices of reference points q\ , q 2 € Q 
and various choices of k G H(qi, q 2 ). 

We are particularly interested in situations when the action of G on Q is 
free, i.e., when the isotropy groups are trivial, H(q) = {e} for any q € Q. Then 
H(qi,q 2 ) are one-element sets, H{q\ 1 q 2 ) — {fc}, i.e., the above fc-element is 
unique. The "labelling" of the right-acting transformation group by elements 
of G depends only on the choice of the reference point. It is clear that for any 
qo € Q, N (H(qo)) — G, i.e., the right-acting transformation group is isomorphic 
with G itself. Choosing some reference point q$ we automatically fix one of these 
isomorphisms. 

The extended affinely rigid body is defined as a system of material points 
constrained in such a way that all affine relationships between constituents re- 
main frozen during any admissible motion. Summarizing the above remarks we 
can formulate a few geometric models of its configuration space. 

1. If we use the standard terms of continuum mechanics based on the affine 
physical and material spaces [M, V, — »), (N, U, — >), then the configuration 
space is given by AfI(iV, M), i.e., the manifold of affine isomorphisms of 
N onto M. Affine groups GAf(M), GAf(A) act on AR(N,M) according 
to the rules ([2]), |3]) and describe respectively spatial and material trans- 
formations (kincmatical symmetries). If we formally admitted singular 
configurations with degenerate dimension, the configuration space would 
be given by Af(A, M), i.e., the set of all affine mappings of N into M, in 
general non-invertible ones. By the way, Af(iV, M) is also an affine space 
with Af(A, V) as the translation space; this is just the special case of (HJ. 

2. If some material origin O € N is fixed, the configuration space may be 
identified with 

Q = QtrX<3 lnt = MxLI(C/,V); 
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the first and second factors refer respectively to translational and internal 
(relative) motion. And again, when singular configurations are admitted, 
LI(Z7, V) is replaced by its linear shell L(U, V). The natural groups of affine 
symmetries act on Q according to (O, {7}, ©, ©• When translational 
motion is neglected, the configuration space reduces to Q[ nt — Ll(U,V), 
or simply to L({7, V) when singular internal configurations are admitted. 

3. When in addition to O G N some linear basis E — (E% , . . . , Ea, ■ ■ ■ , E n ) is 
chosen, i.e., when an affine frames (0, E) is fixed in TV, the configuration 
space becomes identified with 

Q = Qtr x Qtat = M x F(V); 

F(V) denotes as previously the manifold of linear frames in V (n — 
dim V). When we are not interested in translational motion, simply 
Qint = F(V) is used as the configuration space. If, for any reason, singu- 
lar configurations are admitted, we extend F(V) to V n = V X • • • X V (n 
Cartesian factors). 

Just as in the model Q = Mx LI(C7, V) transformation groups act essen- 
tially according to the rules (HJ), (0), ((HI), ([9]), the linear space U being 
replaced by R™ (any choice of E e F(U) identifies U with R n ). More 
precisely, spatial transformations are given by 

A G GAf(M) : M x F{V) 3 (x; . . . , e K , . . .) ~ {A(x); L(A)e K , . . .), 

a G GL(V) : M x F(V) 3 (x; ... ,e K , ...))->■ (x; ... , ae K , . . .). 

The frame (0,E) identifies N with R"; namely, the point a G N is iden- 
tified with its coordinates a K (a) with respect to (O, E): 

Oa = a K E K . 

Therefore, GAf(AT), GL(U) are identified respectively with GL(n,M)x s IR", 
GL(n, K). Their right-hand-side actions on Q — MxF(V) are respectively 
described as follows: 

{13, u) G GL(n,R) x s M" : (x; . . . , e K , . . .) ^ {t eu {x); . . . ,e L l3 L K , . . .), 

(3 G GL(», R) : (x; . . . , e K , . . .) ^ (x; . . . , e L /3 L K, ■ • •)» 

where eu £ V denotes the vector the coordinates of which with respect to 
the basis e coincide with u K , k = 1, n, 

eu = u K eK- 

4. Model with the structure-less material space. This is just the model based 
on orbits and homogeneous spaces, described in some details above. So, 
51 is the (structure-less) set of material points and Inj(f2, M) denotes the 
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set of injections of O into M. The spatial affine group GAf(M) acts on 
Inj(fi,M) through (13]), 

A e GAf (M) : lnj(0, M) 3$k4o$. (14) 

Any orbit of this action may be chosen as the configuration space of affincly 
rigid body. Different orbits are related to each other by non-affine trans- 
formations. More precisely, we usually concentrate on such orbits Q that 
for any $ £ Q C M is not contained in any proper affine subspace 

of M. Therefore, the body is essentially n-dimensional (n = dim M), 
although, obviously, it need not be so in the rigorous topological sense 
(e.g., when f2 is finite the body is topologically zero-dimensional). Let us 
mention, however that there are interesting applications of the model of 
singular affine body. The configuration space Q is then such an orbit of 
(fT4"]l that for any $ g Q the subset $(f2) C M is contained in an affine 
subspace of M of dimension k < n. The right-acting partner of (Tj"4")) is 
then constructed as described above for the general homogeneous space 
(Q, G). Now G — GAf(M) and Q is an orbit of (fT4"|) consisting of injections 
<f> with n-dimensional affine shells of (n = dim M). 

There is another way of fixing the configurations of an extended affine 
body with the structure-less material space. Let us assume that the body 
is non-degenerate. There exists then an (n + l)-element subset B C ft 
of material points such that for any $ £ Q $(£?) C M is not contained 
in any proper affine subspace, i.e., its affine shell coincides with M. Let 
us take the elements of B in some peculiar order (wi, . . . ,uj n +\). Every 
configuration $ S Q is uniquely fixed by position (<j>(ui), . . . , $(w„+i)) of 
the ordered system (u>i, . . . , ui n +i). The current positions $(oj) of all other 
material points u> € Q are uniquely determined by (<&(wi), . . . , $(u n +i)). 
The reason is that all affine relationships between positions of material 
points, i.e., all linear equations satisfied by vectors (ui, u>' be- 

ing arbitrary elements of f2) are invariant during any admissible (affincly 
constrained) motion. In other words, they depend only on to, u>' but are 
independent of $ € Q. 

In this way, configurations are identified with elements of the Cartesian 
product M n+1 = M x ■ • ■ x M ((n + 1) copies of M). When the body 
is non-singular, it is not the total M n+1 that is admitted but its open 
subset consisting of such (n + l)-tuples (y%, . . . , tja, ■ ■ ■ , Un+i), which are 
not contained in any proper affine subspace of M (therefore, the affine 
shell of ua, A = 1, (n + 1)), coincides with M. This means that the 
vectors yTy2, A — 2, (n + 2), are linearly independent. 

For fc-dimensional degenerate body (k < n = dim M) the configuration 
space Q (an orbit of GAf(M)) may be identified with M k+1 = M x • • • x M 
((k + 1) Cartesian factors) or rather with its open subset consisting of 
(k + l)-tuple (yi, . . . , tja, ■ ■ ■ , Uk+i) with fc-dimensional affine shells. 

Non-degenerate ordered (n + l)-tuples y € M n+1 may be interpreted 
as affine bases (affine frames) in M. For any pair of such bases y — 
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(Vi, ■ ■ ■ ,VA, ■ ■ -,y n +i), y' = (y[, ...,y' A ,.. -,y' n+ i) there exists e xactly one 
affine transformation <I> G GAf(M) such that y' A = ^(va), A = 1, (n + 1). 
Any such basis may be interpreted as a reference configuration. 

5. Finally, if some affine frames (0,E), (o, e) are fixed both in N and M, 
these affine spaces become identified with R, and the numerical affine 
group GAf(n,R) ~ GL(n,R) x s R n may be used as the configuration 
space. Spatial and material transformations become then respectively left 
and right regular translations. 

Obviously, all the above models are mutually equivalent and their formal 
utility and practical usefulness depend on the kind of considered problems. For 
example, when describing discrete systems we shall use the model (4) with the 
structure-less material space. In certain microstructural applications and in fun- 
damental physics one uses internal degrees of freedom which are not interpreted 
in terms of composed multi-particle systems and perhaps by some principal 
reasons do not admit such interpretation at all (cf. the concept of spin of el- 
ementary particles). Obviously, then the model (3) based on the manifold of 
linear frames FM is the most adequate one. When one deals with the motion of 
structured bodies in non-Euclidean spaces, this is practically the only adequate 
approach if one wishes to remain within the framework of finite-dimensional 
analytical mechanics [SH EZl EH EE E01 EH EOl ED E2] ■ 

There are some delicate points concerning the connectedness of the configu- 
ration space. Obviously, within the standard continuum treatment the singular 
situations are forbidden, $ and if must be bijections, and e is a linear frame, 
not an arbitrary n-tuple of vectors. Therefore, the genuine configuration space 
is then one of the two connected components of Afl([7, V), LI(J7, V), F(V). Oth- 
erwise one would have to pass through the forbidden " singular configurations" . 
Only the connected components of group unity, i.e., orientation preserving sub- 
groups GAf+(M), GL+(V), GM+(N), GL+(C/), GL+(n,R) are admitted as 
transformation groups. However, it is not the case when one deals with discrete 
affine bodies. Then there is nothing catastrophic in passing through "singu- 
lar" situations when at some instant of time $(^2) is contained in a proper 
affine subspace of M . And there is nothing bad in mirror-reflected configura- 
tions forbidden in continuum mechanics. Obviously, the above description of 
configuration spaces must be modified then, e.g., Q is not any longer the ho- 
mogeneous space of GAf(M); rather, it may be a union of various transitivity 
orbits corresponding to all possible dimensions k < n of affine shells of $(0). 

Analytical formulas will be usually expressed in terms of rectilinear Cartesian 
coordinates x % , a K respectively in M and N. They are fixed by affine frames 
(o,e), (<D,E) ' ^ 

oi~x\x)ei 1 Oa = a K (a)EK- 

Coordinates x l , a K induce parameterization of configurations. Eulerian and 
Lagrangian coordinates (spatial and material variables) are related to each other 
by the formula: 

y l = ^(a) = x l + ^ K a K 1 
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where y l are coordinates of the spatial position of the a-th material point, and 
x % are coordinates of the spatial position <f>(0) of the fixed reference point 
O G N. The quantities [x l , ^k) are labels of $ and may be used as generalized 
coordinates q a , a = 1, n(n + 1) on the configuration space Q. 

The reference point O € N, i.e., origin of Lagrange coordinates (a K (O) = 0) 
was chosen here in a completely arbitrary way. In practical problems the choice 
of O is, as a rule, physically motivated. If there exist additional constraints 
due to which some material point is immovable (i.e., the body is pinned at 
it), then the material origin O usually is chosen just at this point. There is 
no translational motion then. If translations are non-constrained, the centre of 
mass is usually chosen as the material reference point. Let us remember that 
in situations other than continuous medium filling up the whole space, centre 
of mass may happen to be placed "in vacuum". And even if it is not the case, 
the centre of mass is something else than the material point coinciding with it. 

Let the reference mass distribution be described by some positive regular 
measure /i on N; this means that the mass of the sub-body B C N is given by 

H(B) = f d^. 

J B 

Centre of mass C(p) in N is the only point satisfying 

J C(fj,)adfj,(a) = 0; 

the dipole moment of /j, with respect to C(fi) vanishes. Any configuration $ 
gives rise to the ^-transported measure /x$ on M, 

M<fr(£)) - KB), ^{A) - K^HA)). 

The measure /x$ describes the Eulcrian mass distribution on M (current mass 
distribution). The current centre of the mass distribution C(/i$) € M is given 
by the formula: 

J C(n$)ydn$(y) = 0. 

It is well-known that the centre of mass is an invariant of affine transformations; 
by the way, affine transformations may be defined just as those preserving cen- 
tres of mass. Therefore, 

$(<7(/x)) = C(^). 
And besides, for any affine transformation A g GAf(M), 

A(C{p*)) = C(hao$). 

This is not true for non-affine configurations and transformations. 

If the material reference point O is chosen as C(/z), i.e., Lagrangian centre 
of mass, then generalized coordinates (x l ,(p l x) on Q are especially convenient 
because x % are spatial (Eulerian) coordinates of the instantaneous position of 
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the centre of mass in M. The variables ip l x refer to the purely relative (internal) 
motion. 

The physical quantity /i (Lagrangian mass distribution) is at the same time 
an auxiliary geometric object underlying the convenient models Mx 1A(U,V), 
Mx F(V) of the configuration space of the affine body. Practically all over this 
treatment we put the origin of Lagrangian coordinates at the Lagrangian centre 
of mass. 

The centre of mass is defined as such a point with respect to which the dipolc 
moment of the mass distribution vanishes. The monopolc moment equals the 
total mass of the body, 



m = (J-(N) = / dfj,. 

J N 

Higher order multipole moments given an account of inertial properties of ex- 
tended bodies. Thus, the Lagrangian second-order moment J 6 U <8> U is given 

by 

J KL : = f a K a L d f i(a). 



As mentioned, the origin of a^-coordinates is placed at the Lagrangian centre 
of mass C(/i). The above object J KL is algebraically equivalent to the usual 
co-moving tensor of inertia known from the rigid body mechanics. One can 
^-transport it to the physical space M. Obviously, the result J ((f) <E V <g> V is 
non-constant; it depends explicitly on the configuration <3> but only through its 
internal part if, 

J{if)^ = ^ k ^lJ kl . (15) 

It is clear that 

JfoO« = / (f - x l ) M - x>) d^(y). (16) 



Obviously, the tensors J, J(ip) — • J are symmetric and positively definite. 
The most convenient choices of the material reference frames E are those diag- 
onalizing J. 

One can as well define higher-order multipoles, i.e., 



J K i- K < = I a Kl ---a K >dn(a), 



Obviously, in affine motion 

j (if) =<f x Kl ---if i Kl j 

Inertial multipoles of the order I > 2 do not occur in mechanics of affine bodies, 
nevertheless they are useful in other problems of continuum mechanics |13) . 

As yet we have used above only affine concepts, i.e., we remained on the 
ascetic level of Tales geometry. No metric concepts like distances and angular 
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were used. The admissible configurations $ € AfI(iV, M) are homogeneous in 
the sense that for any a e N the placement D a <& e LI([7, V) takes on the same 
value <p = L[<1>]. Nevertheless, it would be incorrect to say that deformations 
are homogeneous, because without metric (Euclidean) geometry there is no de- 
formation concept at all. The only well-defined " deformation" is then violation 
of affine geometry, i.e., non-constancy of the mapping a D a <fr. 

Let us now introduce metrical concepts. The spatial and material metric 
tensors will be denoted respectively by g £ V* <g> V* , r\ € U* ® U* . By definition 
they are symmetric and positively definite. Their contravariant inverses are 
denoted by g e V ®V, rj € U®U, but in analytical expressions we use the same 
kernel symbols; the distinction is indicated by the use of lower- and upper-case 
indices, 

9ij, 9 ij , Vab, V AB ; 9 lk 9kj^5 l v V AC Vcb = S a b . 
For any configuration <f> e GAf(M) we define Green and Cauchy tensors 

G[$] eC/*® U*, G[$] eV*® V\ 

G[$] = <p*g, C[$] = <^- 1 *r ? , 

i.e., analytically: 

G[$] A B = gij^A^B, C[<f>] = VABtp~ 1A t¥~ 1B j- 

In these formulas, as usual, ip denotes the linear part of <&, ip = L[<f>] = D$. For 
general non-affinc configurations these tensors become fields respectively on N, 
M, namely: 



G[$] a = £> a $* ■ g, C{<t>} y = Dy®- 1 * ■ r, 



Analytically: 



dy % dyi da A da B 

= 9H g^A g^B > C V =VAB-^r-g^. 

Obviously, Green and Cauchy tensors are symmetric and positively definite. 
G is built of the spatial metric tensor g and is independent of the material 
metric rj. And conversely, C is independent of g and explicitly depends on rj. 
Therefore, the traditional term "deformations tensors" is rather non-adequate 
here; the deformation concept presumes comparison of two metrics, whereas G 
and C are well-defined even if respectively rj and g are not fixed at all. Let 
us assume they are both fixed and denote the corresponding Euclidean space 
structures by (N,U,—*,rj), (M, V, — >,g). When some configuration $ is fixed, 
any of these two spaces is endowed by two metric-like tensors, respectively, 
G, rj ee U* <g> U*, G, g e V* <g> V*. The Lagrange and Euler deformation tensors 
are respectively given by 

E[*] = \ (G[*]-rj), e[*] = i( 5 -e[$]). 
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They vanish in the non-deformed configurations, i.e., when <!> are isometries. 
Obviously, these are usual (metrically) rigid body configurations. Their man- 
ifold will be denoted by ls(N, 77; M, g) C AK(N, M). The isometry groups 
ls(N,r]) c GAf(A), ls(M,g) c GAf(M) act on Is(N,r);M,g) respectively on 
the right and on the left in the sense of (O, (0]). Obviously, in realistic classical 
mechanics of (metrically) rigid body mirror-reflected coordinates are excluded 
and the genuine configuration space is given by some connected component of 
ls(N, 77; M, g). When orientations are fixed in N, M, this will be the man- 
ifold ls + (N,r];M,g) of orientation-preserving isometries. Physically admissi- 
ble symmetries are given by the connected subgroups Is + (7V, 77) C GAf + (A), 
Is+(M, 5 ) C GAf + (M) of orientation -preserving transformations (obviously, 
these group themselves do not assume any fixed orientation; they preserve sep- 
arately both of them). 

When translational degrees of freedom are neglected, configurations of in- 
ternal (relative) motion are elements of D(C/, 77; V, g), i.e., the manifold of lin- 
ear isometries of {U,rj) onto (V, g). The corresponding spatial and material 
transformations are respectively elements of the subgroup 0(V,g) C GL(V), 
0(t/, 77) C GL(t/); g- and ?7-orthogonal transformation groups. And again in 
classical rigid body mechanics one should restrict ourselves to one of the two 
connected components of 0(U, 77; V, g). When orientations in U, V are fixed, this 
is SO(U, 77; V, g) C LI + ([7, 77; V, g), i.e., the manifold of orientation-preserving lin- 
ear isometries. The connected components are ruled by the proper orthogonal 
groups SO(F, g), SO([7, 77) consisting of isometries with positive, thus plus-one- 
determinants (no orientations in V, U needed for fixing these subgroups). When 
using the matrix representations we describe configurations by elements of the 
orthogonal group 0(ti, R), or in classical problems, by elements of the proper 
rotation group SO(n, R). 

Let us now fix some metric tensors 77, g. The geometric structure becomes 
more rich and certain additional object may be defined. For example, the ma- 
terial reference frame maybe made less arbitrary and more based on physical 
concepts. Let us define two tensor objects built of the inertial tensor J <E U <S>U 
namely J <EU*®U* and J € U(&U* ~ L(C7), analytically given by the formulas: 

7 jCB x b 7a jAC„ 

JACJ — OA , J B = J VCB- 

Obviously, the object J is non-metrical, but J depends explicitly on the material 
metric tensor 77. Now we have a well-defined cigenproblcm in U : 

JE = IE. 

In the generic non-degenerate case there are n mutually distinct eigenvalues I a , 
A = 1,...,ti, and n mutually orthogonal eigendirections. The directions are 
determined by vectors Ea which may be chosen 77-normalized to unity and such 
that the orthonormal frame E = (E±, . . . , Ea, ■ ■ ■ , E n ) is oriented positively 
with respect to the fixed orientation in U . When the eigenvalues I a are ordered 
by convention in increasing order then Ea are unique up to multiplying some 
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of them by minus-one-factors. And the inertial tensor J is then represented as 
follows: 



J= J A E A ®E A . (17) 

A=l 



This is just the best choice of the material reference frame. Obviously, it is no 
longer unique when degeneracy occurs, but these are non-generic situations. The 
extreme degeneracy correspond to the r\ -spherical body, when J is proportional 
to T) 

J AB = m AB . 

Remark: If ip is not an isometry, then obviously the co-moving vectors &a = 
tpEA are not orthonormal in the g-sense, however, they are orthonormal with 
respect to the Cauchy tensor C used as a "metric" in V. And then: 

71 

JM = jAeA ® eA 

A=l 

with the same values J A which occur in (TIT)) . 

Remark: For the sake of economy of symbols we could as well denote J A b 
by J A b, following the convention used in Euclidean and Riemannian geometry. 
But it is not the case with Jab- Writing it as Jab would suggest tjacVbdJ 00 , 
quite incorrectly, this is not the 77- lowering of indices. Typical notational short- 
hands may be misleading. Only for the ?7-manipulation of indices and for the 
contravariant inverse of r\ {r\ AC r\cB — S A b) we can safely use the index manip- 
ulation with non- modified kernel symbols. 

When metric tensors are fixed, we can discuss the problem of interaction 
between rotations and deformations. In any of linear spaces U , V we are given 
two symmetric positively definite tensors: r\ G U* ® U* , G £ U* ® U*, g £ 
V* ®V* , C £ V* ®V* . Obviously, 77, g are fixed whereas G, C are configuration- 
dependent. Just as previously, we can define the byproduct-objects G £ U ®U , 
G £U ®U* ~ L(U), C £V ®V, C £V ®V* ~ L(V). Analytically they are 
given by 

G~ AC s~t ?A pfA ^AC r~i s~iiks~i ?i pti „iks~i 

^CB — OB, B — V l*CB, t> Ujk — j , O j = g Ukj- 

Again the same care must be taken as to the upper- and lower-case indices. 

And now G [</?], C[<p] may be expressed in terms of their rj- and (j-orthonormal 
bases (. . . , F a [tp], . . .), (. . . , f a [tp], . . .) in U, V 



G[<p] = ^2 \ a [ip\F a [<p\ (g) F a [cp], 

a=l 
n 



a— 1 



obviously, (. . . , F a [ip], . . .), (. . . , / a [<^], . . .) are the dual orthonormal bases of 
U*, V*. When there is no danger of misunderstanding, the label ip at A a , F a , 
f a , F a , f a may be omitted. 
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The quantities A a [ip] arc deformation invariants in the sense that they do 
not feel spatial and material linear isometries, 

\ a [ot(p0\ = \ a [tp] 

for any a £ 0(V,g), j3 € 0(U, rj). The more so they are non-sensitive with 
respect to the spatial and material affine isometries (because translations ev- 
idently do not affect them). Unlike this, the Green and Cauchy deformation 
tensors are non-sensitive only with respect to spatial and material isometrics, 

G[ap]^G[p], C[<p0\=C[<p] 

for any a e 0(V,g), f3 E 0(U,rj), but not conversely. F a [tp] e U, f a [tp] G V are 
normalized eigenvectors respectively for G[ip], C[ip}: 

G[<p]F a [<p] = \ a F a [cp], C[<p]f a [<p] = f/„M, 

V(Fa,Fb) = S ab , g{fa,fb) = Sab, 

and they are essentially unique for the non-degenerate spectra. 

The fixed bases (..., F a [tp], ...),(..., f a [tp], ...) , being orthonormal, give rise 
some isometry U[ip] € 0(U,tj;V, g), namely such one that 

U[lf]F a [ip] = f a [ip], 0=1,..., 71. 

The ip itself may be written in the form 

p = U[p]A[p], (18) 

where A[cp] : U — > U is r/-symmetric and positively definite, i.e., 

r\ (A[<^]u, v) = 7] (u, A[tp]v) , 77 u) > 0, 

for arbitrary u,v € U and arbitrary u 5^ in the inequality. Analytically: 

T] AC A[p} c B = r)BcA[ip] c a, 

and the matrix \A[(p] b\ has positive eigenvalues, coinciding, by the way, with 
square roots of deformation invariants A a (which, obviously, are also positive). 
Obviously, (fl8|) is a geometric interpretation of the polar decomposition, and 
analytically, when we put U = V = R n it exactly coincides with the usual polar 
decomposition known from the matrix theory. It may be also alternatively 
written in the form: 

<p = B[<p\U[<p], (19) 
where B[ip] : V — > V is ^-symmetric and positively definite, and obviously: 

B[<p] = U[<p]A[<p]U[<p]- 1 . 
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These are simply the left and right polar decomposition known from the matrix 
theory. It is clear that Green and Cauchy tensor satisfy respectively: 

G[p] = A[f} 2 , C[p] = B[ v ]-\ 

Let us mention, there are also other possible choices of deformation invariants; 
every system of n functionally independent functions of \ a , a = 1, . . . , n may 
be used as a basic system of invariants. Let us remind a few popularly used 
system e.g., 

K a [p] =Tr(GM a ) , o=l,...,n. 

Obviously, 

n 

/C a = ^(A 4 ) a . 

»=1 

Another possibility is the system of coefficients of the characteristic polynomial 
of G[ip], I p [ip]: 



det 



n 

G[<p] A B -\S A B ] = 2(-l) fc /„_ fc MA fc . 



k=0 



Obviously, Iq = 1, and for p = 1, . . . , n, I p is the sum of all possible products 
of p quantities A a with different (but not necessarily disjoint) sets of labels 
a = l,...,n, e.g., 

n 

/i = ^ Ai = Tr (g) , J„ = Ai ■ ■ ■ X n = det (g) . 

i=l 

In the physical three-dimensional case 

^2 = A2A3 + A3A1 + A1A2. 

Let us observe that orthonormal bases (...,/„[<£>],...), (..., F a [(/?],.. .) represent 
formally configurations of two fictitious rigid bodies, respectively in (V, g) and 
(U,ri). They refer respectively to the eigenaxes of the Cauchy and Green de- 
formations tensors, therefore, they tell us how the deformation state is oriented 
with respect to V, U (what are instantaneous positions of deformation ellip- 
soids). Unlike this, deformation invariants contain only the scalar deformation 
about the deformation state (how large are stretchings). The manifold of scalar 
deformation states (parametrized by deformation invariants) may be considered 
as double-coset-space of LI([7, V) with respect to the left and right actions of 
0(V,g), 0(U,r)), 

hrv = 0(V,g)\U(U,V)/0{U,ri), 
or, when reflections are excluded, 

Inv - SO(V, 5)\LI + (^, V)/SO(U, r,). 
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Let us observe that, as usual linear frames F[ip], f[ip] may be naturally 
identified with linear isomorphisms: R[ip] : R™ — > U, L[<p] : K n — > V. As 
they are orthonormal they are linear isometries of (U,t]) onto (R™,£) and of 
(V,g) onto (R™,(5); <5 denotes here the natural Descartes-Kronecker metric of 
R™. Similarly, the dual co- frames F[p], f[<p] may be identified with the inverse 
mappings i?^] -1 : U — > R n , _1 : V — > M"; obviously, they are also linear 
isometries. One can show that <p may be represented as 

<p = L[<p]D[<p]R[<p]- 1 , 

where the linear mapping D[ip] : R™ — > R", i.e., simply a matrix, is diagonal, 

DiM ... 



£>[¥>] = 



... D n [<p] 



= Diag(D 1 [<p],...,D n [<p]). 



If there is no danger of misunderstanding we shall omit the label ip at the 
quantities L, D, R. So, finally, we write 

ip = LDR- 1 

and in this way any internal (relative) configuration p> is formally identified 
with the configuration of two fictitious (metrically) rigid bodies and n purely 
oscillatory degrees of freedom of the stretching state. 

If we put U = V = R™, then the above decomposition (two-polar decompo- 
sition, triple decomposition) is formally obtained from the polar one. Namely, 
for any ip — GL(n, R) one starts from the polar decomposition 

ip = UA, C/eO(n,R), AeSym + (n,R) 

and then A is orthogonally diagonalized, 

A = RDR- 1 , D e Diag(n, R), Re 0(n, R), 

so, finally, 

<p = LDR-\ L = UReO{n,R). 

Unlike the polar decomposition, the two-polar one is non-unique. In the non- 
degenerate case, when all diagonal elements of D are pairwise distinct, this 
non- uniqueness is discrete and controlled by the permutation group inter- 
changing deformation invariants. When degeneracy occurs the non- uniqueness 
is more catastrophic, in a sense continuous, and resembles the singularity of 
spherical coordinates at r = (although, one must say, it is much more compli- 
cated). 

Let us finish with some kinematical concepts. Generalized velocity of an 
affine body is given by the pair (i>,£) G V x L(U,V) consisting of the transna- 
tional velocity v and the internal one £. On a given classical motion R3(n 
(x(t), <p{t)) it is analytically given by the system 

dx l dp 1 a 
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When U = R n , i.e., Q = M X F(V), then velocities are elements of V x V n = 
yn+l_ 

It is convenient to use afhne velocities in the spatial and material (co-moving) 
representations, f2 £ L(V), f2 £ L(J7), namely: 

n = fp = ' = ^ £ = ^ 

They are interrelated as follows: 
Analytically: 



dff^A -lA c>A -1A d f' l B ni i -IB 

^ if j, \l b = <P i ^ ) it j = if A" B<£ j 



Eringen in his micromorphic theory [27( I28[ I29j uses for them the term " gyra- 
tion" . They are Lie-algebraic objects related respectively to the right-invariant 
and left-invariant vector fields and differential forms, 



X[E] v =E<p, X E 



= V>E, 



UJ 



dipip , uj — ip 1 dip. (20) 



In very rough, formal terms we would say that u>, uj are obtained from Q, Q 
via the (^-multiplying. More rigorously, they are respectively L(F)-valued and 
L(U)-va\ued differential one-forms on 1A(U,V). Their evaluations on vectors 
tangent to trajectories just coincide with 17, f2. 

The right and left invariance is meant obviously in the sense of transfor- 
mations (UJ), ([5]). They become right and left regular group translations when 
U = V = R n and Ll(U,V) is identified with GL(n,R). In the above formulas 
for vector fields E and E are respectively fixed elements of L(V), L(Z7). 

Affine velocities are non-holonomic in the sense that there no generalized 
coordinates for which they would be time derivatives. This is due to the non- 
commutativity of the full linear group. 

In continuum mechanics Q may be interpreted in terms of the Euler velocity 
field. Namely, the material point which at a given instant of time passes the 
spatial point y £ M has the velocity: 

v(y) = v + Qxy'; 

analytically, 

v l {y) = v l + n l 3 (y j ~x ] ). 

When the motion is metrically-rigid, i.e., gyroscopic, then affine velocity be- 
comes skew-symmetric with respect to the appropriate metric tensor, 

n l j — —g m gjbfl a = — %*i 

6 A B = ~V AK VBL& L K = -& B A , 
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These objects are the usual angular velocities, respectively in the spatial and 
co-moving representation. In the physical three-dimensional case they are iden- 
tified in a standard way the usual pseudo-vectors of angular velocity; namely, 
in orthonormal coordinates 



9: 



sr 



n i = ip i A n A . 



^ j £ jfe^ 7 

1 

Obviously, e denotes here the totally antisymmetric Ricci symbol and indices 
are trivially shifted with the use of Kronccker-delta (orthonormal coordinates). 

Obviously, in the two-dimensional case (also physically interesting) , angular 
velocities are one-dimensional objects, Q and Q numerically coincide and in 
orthonormal coordinates: 

-1 1 d9 

1 



dt 



-1 

1 



cose 
sin( 



— sine 
cos 8 



For the fictitious rigid bodies corresponding to the polar and two-polar decom- 
positions we also introduce the corresponding angular velocities: 

dU, 



oj = 

X = 

X = 

1? = 

d = 



^u- 1 eso(v, g y cl(v\ 
u-^esoiu^y cl(u), 



dL 
~dt 



L- 1 eSO(V,g)' CL(V), 
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L- 1 ^ g SO(n,K)' C L(n,I 
—R-'eSOiU^ycHU), 
i?- 1 — e SO(n,M)' C L(n,I 



These objects are elements of the of the indicated Lie algebras of orthogonal 
groups. So, ui, x a re g-skew-symmetric, S3, 1? are 77-skew-symmetric, and Xi $ 
are skew-symmetric in the usual Kronccker sense. 

In certain problems it is convenient to use co-moving representation of the 
translational velocity, 



-1A 



dx l 
~dt 



Canonical moments, i.e. linear functional on generalized velocities are pairs 
(p, 7r) G V* x L(V, U), i.e., analytically (. . . ,pi, ,p A i, . . .). Their evalua- 
tions on virtual velocities are given by 



((p, tt), (v, £)) = (p, w> + Tr (ttO = pjV* + p A £ A . 



(22) 
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When U = M" and Q = M X F(V), then canonical momenta are elements of 

V* x V* n = 

The duality between n and £ may be expressed in terms of the duality 
between affine spin and affine velocity. More rigorously affine spin S € L(V) 
and its co-moving representation E £ L(C^) are defined as: 

S = (pK, E = 7T<^. 

Analytically 

One uses also the term "hypermomentum" . The quantities E, E are respectively 
Hamiltonian generators of the transformation groups ^ , §5§ . 

Linear momentum pi generates spatial translations. In many problems it is 
convenient to use its co-moving representation 

p A = Pi<p % A 

which has to do with the material translations. 

One uses also the orbital affine momentum and the total affine momentum 
A, J , given respectively by 

Unlike E the quantities A, J depend on the choice of the origin o G M of affine 
coordinates in M. And J , more precisely Jip) is a Hamiltonian generator (mo- 
mentum mapping) of the centre-affine subgroup GAf(M, o) C GAf(M) (affine 
transformations preserving o) acting through ([2]) . 

The doubled skew-symmetric parts of hypermomenta, 

S l j = E%- - g % gjiE 
L l j = x% -g ik g 3 ix l p kl 
3 l j = L l j + S j, 

are the usual angular momenta: the internal (spin), the orbital, and the total 
ones. They are Hamiltonian generators (momentum mapping) of the corre- 
sponding isometry groups. The quantity 

V A B :=t A B - V AL V B K t K L , 

called by Dyson " vorticity" , is the Hamiltonian generator of the right-acting 
rotation group ([9]). When if is not an isometry, then V is not the co-moving 
representation of S, 

o j T V aV bp j- 

S and V generate respectively spatial and material rotations of internal degrees 
of freedom. 3 and p together are Hamiltonian generators of spatial isometries. 

In a sense, E and E are non-holonomic canonical momenta; non-holonomic, 
because their Poisson brackets do not vanish. 
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The pairing between internal canonical momenta and velocities may be now 
expressed as follows: 



(E, tt) = (E, = Tr (EO) = Tr [YQj = Tr (tt£) . 

Analytically: 

(E,0) = ^E,ft^> = E^fi'j = £ B A d A B - 

Just as f2, n themselves, E, E may be interpreted in terms of right - and 
left-invariant vector fields or differential forms. 

We use the standard conventions of differential geometry according to which 
vector fields with components Z % related to local coordinates z % are identified 
with first-order differential operators 

z = z*° 

Then E, E as systems of vector fields dual to systems of Pfaff forms (|20p (more 
precisely, to L(V)- and L([/)-valued differential one-forms) are given by 

E'j = A~77~j j E A B = ^b^~. (23) 
dipJ A dtp 1 A 

In other words, at the point <p € LI(U,V), the ( fe yi)-th component of E l j 
equals (p l A^ k j, and the ( l c)-th component of E A b equals ip l sS A c- 

Interpreted as invariant forms E- and E-objects become respectively the 
following fields on Ll(U, V): 



Y[F] V = <p- l F, Y 



F 



Ftp' 1 , 



where F, F are arbitrarilly fixed elements of L(V), L(C7), and we remember 
that the dual space L(U,V)* is canonically isomorphic with L(V,U) through 
the formula (|2"2"|) . 

Let us quote the obvious transformation rules of f2, f2, E, E under transfor- 
mations (01, © of internal degrees of freedom: 



a e GL(V) 
P e GL(U) 
a e GL(V) 
£ GL(U) 



fl i— > alia , O i — > 
E i ► aEa -1 , E i— > E 



Now let us quote the basic Poisson brackets. The most important of them, 
namely those involving the above generators, are determined by the structure 
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constans of the linear and affine groups. 

S fe ( } = 5 t iT. k j — S k 3 Y 1 l i, S C d| = 5 C b^ A d - S A d^ C b, 

( rri rrh "1 <ri rrh c/c rri \ ! r, 

\J j, J l) — iJ j — jj i, I h j,2u B f = 0, 

{A'j-.A*,} = 5\K k 3 - 5 k j A i l {t A B ,Vc} = 5 A cPb, 

{j 2 j,Pk} = {A. z j,Pk} = 5 % kPj. 

For any function F depending only on generalized coordinate x % , (p l a we have: 

OF 



{^j,F} = -E i j F = -i P i A 



{£ a b ,f} = -E a b F = -^b 01 



dip 1 a' 



These Poisson brackets are in principle sufficient for obtaining equations of 
motion in the form: 

^ = {F,H}, (24) 



dt 

H donating the Hamilton function. 

It is convenient to introduce in addition to vector fields E l j, E A b and dif- 
ferential one-forms Q A b some others objects, namely, the vector fields H a , 
Ha and differential one- forms 8 a , 8 A , all defined on the configuration space 
Q = MxLl(U, V). In terms of affine coordinates x % , tp l a but result is coordinate- 
independent they are given by 

e a = dx\ H « = ^ e A = <p- 1A idx\ H A = ip i A -£-. (25) 

The following duality relations hold among them: 

/ Q A b,E c d )=6 a d S C b, (0 A b,H c 



A ,E C D )=0, (9 A ,H B )=S A 



B, 



and similarly, 



(Lu\,E c d ) =5 a d 6 c b , (uj a b ,H c ) = 0, 

(e a ,E c d ) = o, (e a ,H b ) = s\. 



Therefore, they are mutually dual fields of non-holonomic frames and co-frames 
on Q. As seen from the structure of Poisson brackets, these fields are geo- 
metrically important. In the theory of principal fibre bundles of linear frames 
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or co-frames they are know as structural fields and standard horizontal fields 

[321 1631- 
Additional important Poisson brackets: 

{ Pa , F} = —H a F, {p A ,F} = -H A F (26) 

for any function F depending only on generalized coordinate. 

Let us finish the above description of classical geometry of degrees of freedom 
with a brief review of symmetry problems underlying the polar and two-polar 
decomposition. They are very important for quantization problems. 

First of all we shall modify slightly our notation. We introduce new gener- 
alized coordinates parameterizing deformation invariants. In many problems it 
is convenient to denote the diagonal elements of D by Q a , a — 1, . . . , n, 

Q a := D aa , X a = (D aa ) 2 = {Q a f . 

And many formulas become remarkably simplified when the logarithmic scale 
is used for parametrizing deformation invariants, q a = lnQ a , thus: 

Q a = D aa = exp (q a ) , X a = exp (2<f ) . 

Deformation parameters q a run over the total real range K. They are fic- 
titious "material points" moving along the real axis. As such they are essen- 
tially identical and indistinguishable; this has to do with the mentioned non- 
uniqueness of the two-polar decomposition. This non-distinguishability is essen- 
tially striking and interesting in the quantized version of the theory. The volume 
extension ratio is given by det D — exp (q 1 + • • • + q n j , thus it may measured 
in a convenient way by the sum (q 1 + • • • + g n ) . It is often convenient to split 
D into the isochoric (incompressible) and the purely dilatational parts, 

D = IA, detA = l, 

The factor / is the linear size extension ratio; obviously, 



I = VdetD = 




The logarithmic measure of this ratio: 

q = \nl = -(q 1 + --- + q n ) 

n 

is simply the " centre of mass" of the mentioned " material points" . The isochoric 
part A depends only on the ratios Q z /Q-* , i.e., logarithmically, on the "relative 
positions" q l — qi . 

The splitting of internal configurations (p into dilatational and isochoric parts 
may be written down as follows: 

tp = W = exp(g)* = exp(g)LAi? _1 , 
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where A as preciously is diagonal and isochoric (det A =1). The term ^ refers 
to the shear-rotational degrees of freedom. It is convenient to use the isochoric 
affine velocities 

dt dt 

They are trace-less, i.e, v e SL(F)', v e SL(f/)' (elements of the Lie algebras of 
SL(V r ), SL(C/)). The total afhne velocities may be expressed as follows: 

= ^ + ^Id v , fUP+^Idu, (27) 

where Idy, Idu are identity transformations in V, U. 

Similarly, the affine spin may be decomposed as follows: 

£ = o-+-Id v , E = ff+-Id v , (28) 
n n 

where a £ SL(V)', a € SL(U)' (traceless) and p is the dilatational canonical 
momentum, 

Tr(E) = Tr (e\ = p. (29) 

This momentum is canonically conjugate to the above q-variable (logarithmic 
size variable). The pairing between velocities and momenta may be expressed 
as follows: 

Tr (SO) = Tr = Tr (au) + pq = Tr (aQ) + pq. 

Poisson brackets for the components of a are based on the structure constants 
of SL(y). The same based for a with the only precise that the signs are reverse. 
The mutual Poisson brackets {a, a} vanish. And obviously, {q,p} — 1, and 
the dilatational phase-space variables have vanishing Poisson brackets with the 
shear-rotational quantities a, a. 

Dilatational canonical momentum p may be interpreted as the total linear 
momentum of the one-dimensional q a -particles. 

The two-polar decomposition identifies (modulo some non-uniqueness) in- 
ternal configuration with the triplets (L; q 1 , . . . , q n ; i?), where (L, R) is the pair 
of rigid bodies (Cauchy and Green deformation tensors principal axes), and the 
fictitious one-dimensional material points q 1 , . . . , q n are deformation invariants 
(in logarithmic scale). The formulas (f2"Tj) suggest us to make use of two possible 
systems of non-holonomic velocities: 

x;...,g a ,-.-;tf) , (x;...,qV ■•;<?)• (30) 

As expected from the (metrically) rigid body mechanics the first subsystem is 
more effective in analysis of dynamical models. 

Similarly, when the polar decompositions are used, we have at disposal the 
following natural systems of non-holonomic velocities: 

(uj,b), (w,A), (uj,b), (31) 
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cf. (HHJ), (JTHJ) . For the qualitative analysis of practically important dynamical 

models the (u), A"j -system is most convenient. 

Non-holonomic canonical momenta conjugate to (|30[) are respectively de- 
noted by 

(P! ■ ■ ■ ,Pa, ■■■;t), (p;... ,p a , . . . ;r) , 

where 

ioGSO(n,M)', T6S0(n,R)', geSO(V,g)', r G SO(U, rj) 1 , (32) 

and p a are canonical momenta conjugate to q a . 

Canonical spin variables of the Green and Cauchy gyroscope and their dual 
angular velocities are considered as elements of the same linear spaces; this 
is due to the natural isomorphisms between orthogonal Lie algebras and their 
duals. So, the corresponding pairings are given by 

((AP,t), = p a q a + ^Tr(px) + ^Tr(rd^ , 

((p,P,T),( X ,ii,#)) = p a q a + ~Tr(p X ) + ^Tr(T&). 
Similarly, the dual objects of (f3"Tj) will be denoted by 

{a, a), (fi,f3), (cr,a), (p, P) , 

where 

a,fie SO(U,r))'; a, fj, G SO(V,g)', 

and a G L(£7), (3 G L(T^) are respectively rj- and ^-symmetric. The correspond- 
ing pairings are given by 

((a, a) , (Q, i) ) = ^Tr (aQ) + Tr (ai) , 

etc., an analogous way for other combinations. 

Let us remind that in the physical three-dimensional case the skew-symmet- 
ric tensors are identified with the axial vectors, e.g., 

X*i = -e'ifcX*, X 1 = -\e % j k X 3 k 

in orthonormal coordinates. For the dual angular momentum quantities we have 
the reversed-sign-convention, e.g. 

~i ^ A k A 

P j = £ jkP , P = ^ e j P k- 

The shift of indices here is meant in the cosmetic Kronecker-delta-sense. Then 
the former formulas are compatible with the standard K 3 - conventions, e.g., 

^Tr(px) = p iX \ 
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and similarly for other angular velocity and angular momentum quantities. 

The quantities p, t coincide respectively with spin S and negative vorticity 
—V. They are Hamiltonian generators of transformations: 

ip^Aip, tp^tpC- 1 , AeSO{V,g), CeSO(£/,r?), 

i.e., in terms of the two-polar decomposition: 

L i— > AL, R^ CR. 

The objects p, f are Hamiltonian generators of transformations 

L^LA, R^RC, i,CeSO(n,R). 

Similarly, if we use the polar decomposition, a coincides with spin S, because 
it generates transformations ip \— ► Wip, W G SO(V,g), i.e., 

U i ^ WU, (p = UA^ WUA = W<p. 

The quantity a is the Hamiltonian generator of 

<p = UA^UKA, KeSO(U,r]). 

Similarly, the objects p generates transformations: 

C E SO(U, rf): <p = BU ^ BUG = <pC. 

Therefore, it coincides with the canonical vorticity V. 
And finally, p generates the transformation group: 

W G SO(V, g) : <p = BU^ BWU. 

Everything said above implies that 

X j — L aX b^ j, 

$ A B = R A a d a b R- lb B , 

w*j = W A Q A B U- 1B 1 

and similarly 

-~a r — lb 



P j — L, aP b^ j, T b — Li a T bK B, 

a j — U ACT B U j, p j — U aP BU j. 



In situations where the variables Q l are more convenient than q % = InQ 1 , 
the canonical momenta Pi conjugate to Q 1 will be used instead pi (conjugates 
of q l ). The relationship is as follows: 

Pi = Pt cxp (-<f) = ^. 
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Having in view applications on the fundamental level, including the atomic 
and molecular structure, we concentrate here on the quantization procedure. 
Because of this, on the classical level we are interest mainly in Hamiltonian 
models. Even if dissipativc phenomena are taken into account, they are consid- 
ered as a correction to the Hamiltonian background. In any case, the primary 
concept is that of the kinetic energy. If we assume that the mechanism of affinc 
constraint is compatible with the d'Alembert principle, then the kinetic energy 
is obtained by restriction of the primary multiparticle kinetic energy to the 
tangent bundle of the constraints manifold. After easy calculations one obtains: 

_ to dx i dx j | 1 d^A dip l B jAB 
-itr + Jint- 2 9ij ^ ^ + 2 9ij ^ ^ , 

where to, J denote, as previously, the total mass and the second-order moment 
of the mass distribution. They characterize respectively the translational and 
internal inertia. It is instructive to quote two alternative formulas: 

T = jG AB v A v B + l -G AB n A K Sl B L J KL , 

T = + \g ij n i k W l J[y] kl . 

Legendre transformation may be written in any of the following equivalent 
forms: 

9T , A 9T . lBA 

Pl = ~£hf = m9ijV ' P 1 = ~dt?A~ = 9lJ " B ' 

PA = -^= mG AB v B , % A B = = G BC n c D J DA , 

BT 

Obviously, it is assumed here that there is no generalized potential depending 
on velocities (e.g., no magnetic forces). Otherwise we would have to replace 
the kinetic energy T by the total Lagrangian. Inverting the above formulas and 
substituting them to the kinetic energy expressions we obtain there formulas for 
the kinetic Hamiltonian: 

T = T tI +7l nt = ^-gVp iPj + ^JABp A iP B j g ij , (33) 

T = J-G ab papb + IjabZ a k X b l G kl , (34) 
2to 2 

T = ^p l p + \j[vh^kVig k \ (35) 

where the "tilda" objects are reciprocal tensors. 

For Lagrangians L = T — V(x, ip), the resulting Hamiltonians have the form 
H = T + V{x,ip). 



31 



The usual kinetic energy quadratic in velocities is geometrically equivalent 
to some Riemann structure on the configuration space, 

1 dn^ da v 

r=-V(g)^^, T = T llv {q)dq»®dq\ 

or, in traditional notation using the arc element: 

da 2 = T^dq^dq". 
For the kinetic Hamiltonian we have: 



T=-T^{q) PtiPv , Y = T^{q) 



2 ^> dq 

Using the previously introduced symbols we can express the metric tensor un- 
derlying our kinetic energy in any of the following equivalent forms: 

T = mg tj dx l ®dx j +g i:j J AB d(p i A®d(p j B, (36) 
T = mG AB e A ®9 B +G AB J KL u A K ®w B L , 
T = mgijO* <g) 9 j + g tj J[(p\ kl Lu i k ® uj J t . 

Similarly, for the inverse metric T we have the following equivalent expressions: 

^ 1 u ® d • . ~ d d 

T = — g 13 — ® — +g lJ J A B- 



m dx l dxi dip' 1 a dpi b 1 

f = -G ab H a ®H b +G ab JklE k a®E l b , 
m 



r 



-g ij Hi « Hj + g ij J[p>} k iE k i ® E l j 



It is clear that the above kinetic energies (metric tensors) are invariant under 
the group of spatial isometries Is(M,g) acting through It is also invariant 
under 0(U,J) acting through i.e., the subgroup of GL(U) preserving J. 
In particular, it is materially isotropic, i.e., invariant under 0(U,r]) when the 
inertial tensor is spherical, J — firj. However, there is no total affine invariance 
either in the spatial or material sense. This kinematical symmetry is broken by 
the tensors <?, J. Therefore, the traditional d'Alembert model of the kinetic en- 
ergy does not belong to the framework of left- or right- (or two-side-) invariant 
geodetic systems on Lie groups or their group spaces, i.e., to the theory devel- 
oped by Hermann and Arnold on the basis of rigid body or incompressible ideal 
fluid dynamics. By the way, with the above metric tensors the geodetic systems 
are non-physical, because they predict the unlimited contraction and expansion 
of the body. And when some extra potential is introduced as a dynamical model 
of deformative vibrations, then, except some very special potential shape, none 
or rather small profit is gained from the group-theoretical model of degrees of 
freedom. 

So, it is a tempting idea to formulate dynamical models unifying two things: 
geodetic description (no potential as far as possible) and affine invariance. It 
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turns out that to some extent this may be successfully done: not only the inertia 
but also interactions are encoded in some afnnely-invariant kinetic energy forms 
(metrics on the configuration space). 

The most general and reasonable class of dynamical geodetic models invari- 
ant under the spatial affine group GAf(M) acting through @ is given by the 
metric tensor: 

T = mr lAB e A ®d B + C b a d c u A b®u c d , (37) 

where C B a D c ar e constant and symmetric in their bi-indices, 

r B d r D B 
>~ A C ~ >~ C A- 

The inverse metric is given by 



where 



T = -v AB H A ®H B + C B a D cE A b <8> E c Dl 

771 



rA K r L C X A s;C 

L, B LJ~ K D = DO B- 



Similarly, the right-invariant metric tensors have the form: 

r = m 9lJ e l <g> e j + K\ l k w*, ® <A (38) 

with similar properties of constants 1Z: 

i k — 7^- k^ i • 



The inverse contravariant metric (underlying the kinetic Hamiltonian) is given 

by 

f = -g^Hi ® Hj + IVhE 1 , g> E l k , WhK k x a h = 5 j b S a i. 

m 

The corresponding explicit expressions for kinetic energies and kinetic Hamilto- 
nians are: 

(39) 



(40) 



In certain problems it is convenient to use another equivalent expressions for 
the translational parts. They are respectively given by 

777 I — ■ 

T tr = T tI = —Cv PiPj , (41) 
2 2m 

T tr = ^G AB v A v B , T tr = ±G AB p A PB. (42) 

2 2m 



T = 


777 ^A-R 

—Vabv a v b 


i ^ rB D c,A oC 
+ ^ A C ^ L D ' 


T = 


* „AB^~ i 1 r B D ?A ?C 

- — 7? paPb + t:>~ a b^ d 




2m 


2 




m i i , 




T = 


—g tj v v 1 + 


i di 








T = 




- —!<■ i c?-> d- 
2 
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Let us observe that expressions for GAf(M)-invariant may be interpreted in the 
following way. No metric tensor g in the physical space is assumed and even if it 
exists (as it does in reality) it does not enter the kinetic energy expression (if it 
did, afRne symmetry would be broken and restricted to isometric one). The role 
of the metric tensor in contraction of tensorial indices is played by the Cauchy 
tensor C. 

There is no kinetic energy model invariant simultaneously under spatial and 
material affine transformations. More precisely, any symmetric twice covariant 
tensor field on Q = M X Ll(U, V) must be degenerate. This is due to the 
very malicious non-semisimplicity of the affine group. However, if we neglect 
the translational motion, then there exist internal metric on LI([/, V) invariant 
under both spatial and material (homogeneous) affine transformations (J4j) , d5j) . 
They are given by 

rf nt = Au K L ® Q L K + Buj k k <8> Q L L (43) 
= Auj k i <£> J k + Buj k k ® uj 1 i, 

A, B denoting constants. Their inverses have the form: 



1 „i „, B 

k 



^-E k , ® E l k - — — E k k ® E l , 



A 1 K A(A + nB) ' 

Such a metric is never positively-definite. The reason is that SL(n,R) is non- 
compact and semisimple. rP nt becomes the usual Killing metric when A = 2n, 
B = —2. But this is the pathological situation, because r? nt is degenerate for 
A/B = —n (due to the dilatational centre in SO(n,R)). 
The corresponding kinetic energies are given by 

T° t = ^Trf^W^fTrtt) 2 (44) 



2 V J 2 
= ^Tr(fi 2 ) + |(Trfi) 2 , 

^ - i Tr ( £2 )-2l(^)( Tr£ ) 2 ^ 
- ^^ ^AiAlnB) ^ ' 

(46) 

The B-controlled term in T mt above is a merely correction. The main term 
(A-controlled one) has the hyperbolic signature (n(n + l)/2 +,n(n— l)/2 — ), 
where the "plus" contribution corresponds to the non-compact dimensions and 
the "minus" one to the compact dimensions in GL(V), GL(U). This is the 
highest possible symmetry of Ti nt , an affine counterpart of the spherical top. 
One is rather reluctant to non-positive "kinetic energies". However, one can 
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show that in the above model the lock of positive definiteness is not essentially 
embarrassing; on the contrary, the negative contributions may encode the at- 
tractive part of the deformation dynamics. By the way, the same effect may be 
obtained within the framework of positive Riemannian structures on Q, when 
we use a slightly modified version of pi)) . 

Let us observe that translational kinetic energy (|3"5|) , (piTj) is affinely- invariant 
in the physical space and isometry-invariant in the material space. And con- 
versely, (HO)) . is isometry-invariant (homogeneous and isotropic) in the phys- 
ical space and affinely invariant in the material space. This focuses our atten- 
tion on Riemannian structures on Q — M x LI(J7, V) invariant under the spatial 
affine group GAf(M) and the group of material isometries Is(iV, rf); the opposite 
models are those invariant under spatial isometries Is(M, g) and material affine 
transformations GAf(iV). The corresponding metric tensors are respectively 
given by 

T = mr} KL P< ®^ + Ir )K LV MN 3 K M®u L N + rl lV 
T = mg^^O* +Ig ik gl l U i j ®U k l +Tl A: 

where the constants /, A, B are generalized moments of inertia. 
The corresponding contravariant inverses have the form 

f = ±r ) KL H K ®H L + lr lKL r 1 MN E K M®E L N 
m J 

+ Le k l ®E l k + Le k K ®E L L , 
A B 



r = -g i iH i ®H j + l,g ik gi l E i j ®E k l 

771 I 
1 „ 1 



t- -^E l j®E^ + ^E k k ®E l u 
A B 

where the inertial constants /, A, B are given by 

I=j{l 2 -A 2 ), A=i(A 2 -/ 2 ), B = -±(I + A)(I + A + nB). 
The corresponding kinetic energies are explicitly given by 

T = T tr + T int = ^ mB v A v B + I - mL n K M ^ L NV MN (47) 
A_ / Ao \ B / ^x 2 



T = T tI + T int = jgijvV + ^feQ^OV' (48) 

+ ^Tr (n 2 ) + | (Tr f!) 2 . 

Obviously, the last two terms in both expressions coincide because Tr (J7 P ) = 
Tr (p' P ^j for any natural p. 
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The corresponding kinetic Hamiltonians have the following form: 

T = T tr +T int = ^- 11 AB p A p B + L mL t K M ^ L N7 1 MN (49) 
Zra 21 



Ltt (e 2 ) + X (Tr e\ 
■A V / or V / 



2A V ' 2B 
T = T tI + T int = -LgVpiPj + ^Lg^jX^g* (50) 

+ ^Tr (S 2 ) + -L (Tr E) 2 . 
2A V ; 2B 

Let us observe that the metrical (g- and 77-dependent) parts of kinetic energies 
may be alternatively written down in terms of the Cauchy and Green tensors: 

^-C^ViVi + -iCtfE^EJ'jC* 1 , -^—G^^paPb + Lg ab £ a c V b d G cd . 
2m 21 2to 2/ 

It is important that in a certain open rang of triples (/, A, £?) G M 3 the above 
kinetic energies are positively definite and at the same time they have all ge- 
ometrical and analytical advantages of invariant geodetic system on the group 
manifolds. 

One can show that the spatially affine and materially metrical model (|4"9"]) . 
or more precisely, its internal part, may be expressed as follows: 

1 /--A 1 / ~\ 2 1 „ „2 



1 I , , 



where 

a = I + A, = -~(I + A)(I + A + nB), (i = j (I 2 - A 2 ) , (52) 
and ||y|| denotes the magnitude of the vorticity, 

\\V\\ 2 = -^Tr (V 2 ) . 
Denoting the k-th order Casimir invariant built of generators by C(k), 

C(k) = Tr = Tr (£ fe ) , 

we can write simply: 

T int = ^-C(2) + ±-C(l) 2 + ±\\V\\ 2 . (53) 
Za Zp Z\x 

Similarly, for the spatially metrical and materially affine model (|50p we have 
^^C(2) + A C (l) 2 + i-||,|| 2 , (54) 



3G 



with the same as previously convention concerning the magnitude of spin: 

\\S\\ 2 = ~Tr(S 2 ). 

Let us note that ||V||, | ] -S' 1 1 are simply second-order Casimir built of vorticity 
and spin. 

For the model (|43[) afhncly invariant both in the physical and material space 
we have: 

C = ^cm + iM(n l^ cw. (55, 

It is very convenient to separate dilatational and incompressible motions, espe- 
cially when afhnely invariant kinetic energies (metrics on Q) are used. 
One can easily show that for the afhne-afhne model (H3")l we have 

TL = ^Tr (, 2 ) + n ( A + nB ) e = T n h + r O u (56) 

cf. (|27p . Performing the Legendre transformation we obtain 

^. = ^> 2 > + MaW 2 = T - S + 7;S ' (57) 

cf. (|28|) . (|29|) . Similarly, for the affine- metrical and metrical-affine models one 
obtains respectively the following expression: 

T " = W^A) C ^ {2) + 2n(I + 1 A + n B / + 2^TW) W ^ ' (58) 

Tint = WTA) Csl ^ (2) + 2n (I + A + ub/ + 2(P- W) ^ ' (59) 
where CgL(„)(A:) are Casimir invariants built of a, 

C SL(n) (fc) = Tr( ( r fc ) =Tr(a fc ). 

Let us observe, the only difference is that concerning the last, third term. 
And the both expressions reduce to ([57} when we substitute 7 = 0. And con- 
versely, they may be obtained from (|57[) be replacing: An (I + A) and intro- 
ducing the mentioned terms, 

^aff-mct = j-0 t[A _ j + A] + 2 *_ \\V\\ 2 , 

^ct-aff = „ j + A] + \\ S f . 

Our philosophy is to base the dynamics as for as possible on geodetic affinely- 
invariant models. In particular, geodetic afhne-isometric and isomctric-affine 
models are of special interest. They are "as affine as possible" and at the 
same time compatible with the positive definiteness demand. Nevertheless some 
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models with potentials are still of interest, and, for non-affine models they are 
just unavoidable. So, we shall consider also potential models 

H = T + V, 

where V depends only on the configuration variables (x,cp). What concerns 
inertial properties we concentrated on highly-symmetric models; in any case 
they are always spatialy- and usually materially-isotropic (one can be general 
in formulation, but no so much in effective analysis). It is natural to assume 
that the potential energy V is compatible with these invariance properties of 
the kinetic term. So, V is invariant under internal spatial rotations if and only 
if it depends on ip through the Green tensor G. It is invariant under material 
spatial rotations if and only if it depends on ip through the Cauchy tensor. 

And finally, V is both spatially and materially isotropic in internal degrees 
of freedom if and only if it depends on (p only through the deformation in- 
variants, parameterized, e.g., by q , . . . ,q n . There is a very important special 
case when V is invariant under the volume-preserving groups SL(V), SL(U). 
This means that it depends on ip through the determinant det ip. If we use 
the logarithmic scala of deformation invariants, this means that V is function 
of q = (q 1 + ■ ■ ■ + q n ) /n, the "centre of mass" of logarithmic deformation in- 
variants q l , i = 1, . . . , n. In kinetic energy models {57)) , l[58|). (|55j) dilatational 
and shear-rotational degrees of freedom (incompressible motion) are mutually 
orthogonal; there is no interaction between them. This suggests us to concen- 
trate also on adapted potentials where these degrees of freedom are explicitly 
separated, 

V(q 1 ,...,q n )=V da (q) + V sil (...,q i -qi,...); (60) 

the labels " dil" and " sh" refer respectively to " dilatation" and " shear" . The 
most natural scheme for 14h is that of " binary interactions" between deformation 
invariants 

V* <,')■ (61) 

For isotropic models H lnt = % n t + V (q 1 , . . . , q n ) with T int given by ([5T]) , 
(|58|) the vorticity V is a constant of motion and the third term in (|58|) has also 
the vanishing Poisson brackets with all terms of (|58"|) . The structure of Poisson 
brackets and equations of motion (|24[) implies that the evolution of variables 
q a , ruled by the above Hamiltonian Hi nt , is the same as one ruled by 

Hl t =T^[A^I + A]+V(q\...,q n ), 

where T^ t [A h- > / + A] is obtained from (|45|) T® lt by substituting (I + A) instead 
A. The difference occurs only in degrees of freedom ruled by SO(V, g), SO(U, rj), 
i.e., in the time evolution of quantities L, R describing the orientation of prin- 
cipal axes of deformation tensors C, G. If V depends only on the dilatational 
invariant q, than the total motion in Q is a direct product of two independent 
things: the geodetic incompressible motion and the autonomous dynamics of 
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the g-variable. The deviator 



a 



is then a constant of motion. The general solution for geodetic models based 
on T^ t (|45j) is explicitly given by exponential mapping. Roughly speaking, it is 
produced from initial conditions by one-parameter subgroups of GL(V), GL(J7). 
And it may be shown on the basis of the properties of matrix exponents that 
for the incompressible geodetic affine-affine model 



(with constrains q — 0) the general solution contains an open-subset of bounded 
(oscillating) motions and an open subset of unbounded (escaping, dissociated) 
motions. When dilatations are allowed, then for any Hamiltonian 



with 7J° t given by (|57|) and V(q) stabilizing dilatations, there exists also an open 
subfamily of bounded motions (and an open subfamily of unbounded motions 
if supV^ < oo). The same remains true for the general geodetic affine- metrical 
model ([47)) , |49|) . (|53|) with incompressibility constraints q = 0, and similarly, 
without such constraints but with dilatations - stabilizing potential V(q). 

The same arguments may be applied to dilatationnally stabilized geodetic 
models in (j48|) . (f50|) . (|54|) invariant under 0(V 7 g) x GL(U) or purely geodetic 
isochoric models with the symmetry group 0(V, g) x SL(£7) (materially special- 
affine and spatially material models). On the level of state variables T, A b, q a 
the time evolution is exactly identical with that based on the affine-affine model 
of 7i nt (again with A in ([50]) replaced by I + A). 

Let us summarize the main message. Incompressible amnc-affinc, affinc- 
material and metrical-amnc models may encode the dynamics of elastic vibra- 
tions without any extra potential used, because their general solutions contain 
open subset of bounded motions. When no incompressibility constraints are 
imposed, the same may be achieved by introducing some dilatations-stabilizing 
potential V(q), e.g., some potential well oscillator V(q) — (k/2)q 2 , etc. 

The bounded or unbounded character of motion has to do only with the 
time evolution of g a -variables, and from this point of view the mentioned three 
models are essentially identical. The difference appears only on to level of L, R - 
degrees of freedom, but these gyroscopic variables with compact topology cannot 
influence the property of trajectories to be bounded or escaping. 

To finish this classical description we describe everything in terms of the two- 
polar decomposition. It is convenient to combine the non-holonomic canonical 
momenta p, r in the following way: 




H = 7T° t + V(q) 



M := -p- t, 



N := p- ?. 



(62) 
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This provides a partial diagonalization of the kinetic energy and elimination 
of certain interference terms. Namely, after some calculations one obtain the 
following expressions for Casimirs: 

Cm-Tv 2 I 1 V {Mab)2 1 V {Nab)2 (63) 



and, obviously, 



1fi ^ dh 2 qa ~ qb 16^rh 2 ^ 

a,6 b 2 a,b Ln 2 



C(l) =p = £>«- (64) 



The first term in C(2) may be suggestively decomposed into the "relative" 
and the over-all (" centre-of-mass" ) parts: 

T7- (Pa - Pb) 2 + — . 

a, 6 

This enables one to separate the incompressible and purely dilatational parts. 

It is interesting that C(2) and therefore the kinetic energy itself has a char- 
acteristic lattice structure known from the theory of one-dimensional many- 
body system. Here the logarithmic deformation invariants q a are positions of n 
indistinguishable fictitious "material points". Unlike in the usual Sutherland, 
hyperbolic-Sutherland and Calogero-Moser lattices where all binary coupling pa- 
rameters were identical, now the quantities M l j, N l j are not only non-identical, 
but also non-constant. Moreover, they are state variables subject, together with 
other ones, to some closed system of evolution equations (|2~4")) . where F runs over 
the quantities q l , Pi, L, R, M l j, N' l j. 

Obviously, one should substitute to (|24|) the following basic Poisson brackets 



U a , Pb } = s\, 

{q\ M c d } = { Pa ,M c d } = {q a ,N c d } = { Pai N c d } = 0, 
{M ab ,M cd } = {N ab , Nd} = McbQai - M a d9cb + M ac 9db - M db g ac , 
{M ab ,N cd } = N cb g ad ~ N ad g cb + N ac g db - N db g aCl 

where the shift of indices is meant in the g-sense. Obviously, usually Cartesian 
orthonormal coordinates are used and then simply g ab — S ab . The Poisson 
brackets for M = —p — r, N = p—r follow from the following ones for p, r: 

{Pab, Pcd} = —Pcbgad + Padg c b — Rac9db + Pdfcffac, 

{T abl T c d} = —T cb g ad + T ad g cb — T ac g db + T db g ac , 

{Pab,Tcd} = 0. 

And these brackets are based on the structure constants of SO(n, R), because p, 
t (similarly like S, —V) are corresponding Hamiltonian generators of SO(n, R). 

In affine kinetic energies (|57p . ([55]) . ([55]) the second-order SL(n, R)-Casimir 
invariant has the form: 

CsHn) = ^(Pa~Pb) + ^LtW - T^-^F^- ^ 

a,b a,b bn 2 a,b U1 2 
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It is seen from (J53J|, (O, J55J), JSTJ), (5SJ), |J5SJ), (JBSJ, (53]), ((Ml) that the M-term 
describes some effective centrifugal repulsion of deformation invariants, whereas 
the TV-term is a model of " centrifugal attraction" between q a - n particles" . This 
has nothing to do with any potential V (q 1 , . . . , q n ); this attraction is due only 
to the negative contribution to the affine-affine kinetic energy (|F7)) . In this way 
an apparently "embarrassing" turns out to be just desirable. 

The affine- metrical and metrical- affine models (l58l) . (|59| may be respectively 
written as follows: 

^nf- m0t = 7 £~' S+ 2(P-A*) l|F " 2 ' (66) 

/rmet-aff' _ /7-aff— aff ( 11 oil 2 //^T^ 

J int - J int + 2 (72 _ ^2) 11*11 > 

where 7T^f ~ aff here is just (JS7J) but with A replaced by (I + A). Let us repeat 
the explicit formula for 7T^ t ~ a in terms of the two-polar parametrization and 
the shear-dilatation splitting: 

T-afl-afF 1 C \ 2 1 1 (M ab ) 

T ** = Z(7Ta)^ ( } + WTTa) ^ (68) 



1 {N ab f , p 2 



32(7 + ^1) ^ c h 2 £lzV 2n(/ + A + nS)' 

It is interesting that for Hamiltonians of the form 

H = T+V(q\...,q n ), 

with potentials depending on deformation invariants only, the models (|66[) . (|67p . 
(|68|) give exactly the same evolution equations for the system of state variables: 

;....</"....:..../<„....: \/", : V", ). 

This follows from the basic Poisson brackets quoted above. The only distinc- 
tion between these three model appear on the level of variables L, R, i.e., the 
principal axes of the Cauchy and Green deformation tensors. These degrees of 
freedom have compact topology, thus they do not influence anything concerning 
the bounded or unbounded character of motion. 

It is roughly seen from ([66]) . ([67]) . ([68]) and may be rigorously shown that; the 
incompressible sector of our state variables admits an open family of bounded 
motions and an open family of unbounded ones even in the purely geodetic 
models (without potential). And this is interesting because invariant geode- 
tic systems on Lie groups (SL(n, R) this time) may be successfully analyzed in 
terms of the exponential mapping and special functions on groups. Obviously, 
the dilatational sector violates these nice features. Without potential energy 
the dilatational parameter q moves uniformly in time and the total motion is 
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unbounded. The only bounded solutions q = const are exponentially unsta- 
ble on the level of physical ip- variables. Therefore, the "maximally geodetic" 
affinely-invariant systems have the form: 

H = T + V(q), 

where T stabilizes dilatations. There is no interaction between dilatational and 
shear-rotational degrees of freedom. Dilatational parameter q is subject to the 
one-dimensional dynamics ruled by the Hamiltonian 

2n{l + A + no) 

The same is true in a more general situation when the potential energy de- 
pends also on the shear variables (non-geodetic models) and has the explicitly 
separated form 

V(q\ ...,q n )= V dil (q) + V sh (...,q i - q\ . . .) 

cf (|60|) . (f6Tj) : usually the effective models of V s h will have the binary structure 
(Hill- 
Finally, let us quote the two-polar representation of the doubly-isotropic 
d'Alembert model: 

T. L>T P 2 , 1 V- (Mab) 2 , 1 \- {N ab f 

int ~^Y a w - Q b f 81 h (Q a - Q b f ■ m 

As mentioned earlier, on the purely geodetic level it would be completely 
non-physical. Here it is seen explicitly that 7i nt is purely repulsive on the level 
of Q-variables. All realistic models should be based on some potential term, 

H = T int + V{Q\...,Q n ) 

The binary structure of % nt resembles the Calogero-Moser lattices. And in fact, 
the general scattering solution of the Calogero-Moser chain is a subfamily of the 
general solution of the geodetic model (|69|) . 



2 General ideas of quantization 

After all above classical preliminaries we can formulate the general ideas of 
quantization. We practically restrict ourselves to models based on Riemannian 
structures in the configuration space and on kinetic energies quadratic in ve- 
locities. Let us mention in connection with this the idea of Capriz [S] about 
kinetic energies of more general type, i.e., non-quadratic ones. Such models in 
fact appear in relativistic problems and may be useful in complicated problems 
of condensed matter theory, defects dynamics, etc. However, it may be very 
difficult to use them in quantization problems, because they may need the use 
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of pseudo-differential operators; this may be hopelessly difficult in curved con- 
figuration spaces. So, we remain within the traditional Schrodingcr framework. 

Let us assume that the classical problem is based (as above) on the Riemann 
structure T, i.e., on the kinetic energy form 

r-lr ( ^<W_ 
2 <* vW dt dt ' 

or, if canonical language is used, on the kinetic Hamiltonian 

T=V<^) Pa jv 

The Riemannian volume element is given by 

d^r(q) = yj\det [IV] IV • • • dq f . 

Quantum-mechanical formulation is based on the Hilbert space L 2 (Q, /Ur) of 
C-valued functions with the scalar product 

(*i|* 2 ) = J ^i(q)^ 2 (q)d^r(q)- 
Quantum operator of the kinetic energy is given by 

h 2 

T = --A(T), 

where h is the " crossed" Planck constant and A (r) denotes the Laplace-Bel- 
trami operator of T: 

a (r) = -j= ]T d » (VFlr^dJ) = r""v„v w . (70) 

Obviously, V denotes the Levi-Civita covariant derivative induced by the T- 
metric. This means that the quantum kinetic energy is obtained from the clas- 
sical one by the formal replacing of p M in T by the operator p M = (h/i)V^. 
Parallel transports preserve T and i/jrj, thus, p p is formally self-adjoint in 
L 2 (<2,^r)- When the classical problem is non-geodetic and based on some po- 
tential V(q 1 , . . . , qf), i.e., on the Hamiltonian H = T+ V, then the correspond- 
ing quantum Hamiltonian is given by 

H = T + V, 

where V denotes the operator multiplying wave functions by the potential V, 
i.e., = usually we do not distinguish them graphically. Velocity- 

dependent generalized (magnetic) potentials are not considered here. 

Strictly speaking, from the very principal point of view wave functions are 
not scalars but scalar densities of weight 1/2, and the squared moduli 5 n J' are 
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scalar densities of weight one. But in all realistic models, and we do not go 
outside this scope, some Riemann structure is used and all tensor densities 
are factorized into tensors and standard densities built of T. In particular, 1/2- 
densities "J describing pure quantum states are factorizing as ^ = \/\T\, where 
\I/ are just the aforementioned scalar wave functions. 

Let us also mention that despite some current views, the one-valuedness of 
wave functions is nota fundamental assumption of quantum mechanics. There 
are at least some situations where multivalued amplitudes seem to be acceptable. 
First of all, it is so when the configuration space Q is multiply-connected and has 
a finite homotopy group. Then it is natural to define the wave functions on the 
covering manifold Q. They need not project onto Q as one-valued amplitudes 
but it seems natural to demand that, according to the statistical interpretation, 
the squared moduli \&\I> are uniquely projectable. This has to do with the 
projective representation. And just such situations are interesting in our model, 
where the configuration spaces of rigid and affmcly-rigid bodies in dimensions 
n > 3 have two-element homotopy groups. This point was stressed, e.g., in 
[H [31 [H O 331 HS] , where the possibility of doubly- valued wave functions for 
quantized rigid body was pointed out. 

Literally performed calculations of Laplace-Beltrami operators are usually 
very difficult and the result is rather non-readable. It is much more convenient 
to use directly the operators H^, Hi, E A B , E l j introduced formerly and the 
classical expressions for kinetic Hamiltonians based on these quantities. Then 
we can define easily quantum operators representing the corresponding physical 
quantities, e.g., p;, p A i, 

Hilbert spaces may be constructed without calculating the complicate coor- 
dinate expressions for the metric tensors and their densities \/|Tf. Namely, 
our configuration spaces may be in a sense identified with Lie groups (more 
precisely their group spaces), therefore we can simply use Haar measures, which 
are explicitly known and given by simple expressions. We are usually dealing 
with left- and right-invariant metrics F, thus, the corresponding measures /xp 
are also invariant, and just coincide with the invariant Haar measures, because 
the latter ones are unique (modulo normalization). 

First of all, let us observe that our configuration space Q = Mx Ll(U, V) 
as an affine space (with the translation space Vx L(U, V)) is endowed with the 
natural Lebesgue measure a unique up to normalization. Fixing metric tensors 
g G V*®V* , rj € U*®U* and some adapted Cartesian coordinates (orthonormal 
with respect to these tensors) x l , a K , (f l K, we can normalize I as follows: 

da(x, ip) — dx l • ■ ■ dx n d(p l i ■ ■ ■ d(p n n . 

When translational degrees of freedom are neglected, we use the usual Lebesgue 
measure I on LI(Z7, V) as an open subset of the linear space L(£7, V): 

dl((p) = dip 1 1 ■ ■ ■ dip n n . 

These measures are invariant under translations in the affine space Mx LI(U, V) 
and under spatial and material isometries. They are, however, non-invariant un- 
der spatial and material affine transformations. To achieve the affine invariance 
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we must use the following Haar measures a, X on Q = Mx LI(£/, V) and Q m t = 
U(U,V) induced from the affine group GAf(n,R) ~ GL(n,R) x s R" and the 
linear group GL(n, R): 

da(ip, x) = (det (fi)~ n ~ 1 da(x, ip) = (det (p)~ n ~ 1 dx 1 ■ ■ ■dx n dip 1 1 ■ ■ -d(p n n , 
d\((f) = (det <p)- n cU(<p) = (det(p)- n d<p 1 1 • 

Expressing the measure I in terms of the two-polar decomposition tp = LDR^ 1 
we obtain 

d\((f) = dX(L; q a ; R) = ]J |sh (q l - q j ) \ dq 1 ■ ■ ■ dq n d^{L)d^{R), 

where fi is the left- and right-invariant Haar measure on the manifolds of linear 
isometries Us(R n ,6;V,g), Us(R n ,S;U,ri). Obviously, when LI(U,V) is identi- 
fied with GL(n, R) and the mentioned manifolds of isometries are identified with 
the orthogonal group SO(n, R), then \i becomes simply the literally understood 
Haar measure on SO(n, R). As the manifolds Lis and the group SO(n, R) are 
compact, the measure /i may, although need not, be normalized to unity (the 
manifold volume equals the unity). In certain formulas it is convenient to use 
the symbol 

P A :=ni sh (^-9 j )l' ( 71 ) 

thus, 

dX(ip) = Pxdq 1 ■ ■ ■ dq n dfj,(L)dfj,(R). 

One can also obtain the following convenient expression for the Lebesgue mea- 
sure I: 

dl = PidQ 1 ■ ■ ■ dQ n d^(L)dn(R), 

where 

p i = n (q 42 - q' 2 ) = n (q* + - ® j ) > ( ?2 ) 

and, as we remember, 

Q a = cxp (q a ) . 

The Haar measure on the internal configuration space of the isochoric (incom- 
pressible) affinely-rigid body may be expressed in terms of the Dirac distribution 

dX SL (p) = P x (q\ . . . , q n ) 6 (q 1 + ■ ■ ■ + q n ) dq 1 ■ ■ ■ dq n dfi(L)dfi(R). 

Our quantum- mechanical models will be based on Hilbert spaces L 2 (Q,a), 
L 2 (Qint)0' L 2 (Q,a), and L 2 (<3i nt ,A). Obviously, for afnnely-invariant models 
L 2 (Q,a), L 2 (Qint, X) are more convenient. Similarly, for the usual d'Alembert 
models L 2 (Q,a), L 2 (Qi nt ,Z) are more natural. Nevertheless, it is a matter of 
convenience; one should stress that both types of models may be formulated in 
terms of any of these Hilbert spaces. 
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The spatial and material actions of GAf(M) and GAf(iV) ©, © on the 
configuration space Q preserve the Haar measure a. Similarly, , preserve 
the Haar measure A on the internal configuration space Qi n t- On the other hand, 
except isometries, they do not preserve the usual Lebesgue measures on affine 
spaces, i.e., a, I. The latter ones are invariant, however, under the usual affine 
translations given analytically by 

(. . . , x\ , tp* A , ...) h+ (. . . , x l + C <p' A + Ca, ■ ■ ■) 

just the usual additive translations in Mx L(U, V) and L(U, V) as affine spaces. 
(Remark: on Q and Q- m t, when L(U, V) is replaced by its open subset LI(f, V), 
then these translations act only locally.) On the other hand, these additive 
translations in general do not preserve the Haar measures a, A. 

All the mentioned groups act argument-wise on wave functions. When they 
preserve the measure on Q or Qi u t, the resulting transformations of wave func- 
tions preserve the corresponding L 2 -spaces and are unitary, i.e., they preserve 
the scalar products too. 

Let us quote explicitly some expressions, at least to fix the notation used 
later on. 

For any A £ GAf(M) we define the operation A which transforms the wave 
function W : Afl(N, M) — > C into such one that 

(A$) ($)=*(Ao$). (73) 

Similarly, for any A E GAf(A) we define the operator B such that 

(Bf)($) = $($oB). (74) 

If translational degrees of freedom are neglected and we deal with wave functions 
* : LI(U, V) -> C, then for any A G GL(V), B G GL(U) we define 

(A*)( V ) - V(A<p), (75) 
(B*)( V ) = 9(<pB). (76) 

Obviously, for any A, B the operators A, B are unitary in L 2 (Q, a), L 2 (Q int , A), 
because the measures a, A are invariant under regular translations. Unlike this, 
they are not unitary in L 2 (Q, a), L 2 (Qi nt , A), unless A, B are volume-preserving 
mappings, i.e., elements of SAf(Af), SAf(iV), SL(V), SL(U) (more precisely, 
unimodularity is sufficient, i.e., detL(A) = detL(B) — ±1). Obviously, the dif- 
ferential operators (vector fields) H a , E l j defined in (|25|) . ([23| are generators 
of unitary groups defined in (f73|) . ([75]) . therefore, they formally anti-self-adjoint 
in L 2 (Q,a), L 2 (Qi n t,A), or rather in the subspaces of smooth functions. Be- 
ing non-bounded (non-continuous) they are not anti-Hermitian in the rigorous 
mathematical sense. However, they are so in rough terms used in physics. They 
possess anti-Hermitian extensions. The following differential operators: 

H h d h h a d 

Pa = -H a = -3—, S b = -E b = ~<P k » h 

l l OX l l Oif • K 
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are formally Hermitian. They are respectively quantum linear momentum and 
hyperspin operators. One can also introduce the operator of the total affinc 
momentum (hypermomentum) 

3\ = x a P a + S a fc = A a b + £%. 

The ordering of non-commuting operators meant just as written above; it follows 
from their geometric nature as group generators. Obviously, the coordinate 
operators are defined in the usual way, 

(x°*) (x, If) = X a ^(x, <f), (<p°JC*) (X, If) = f a K^(x, tp). 

If we define quantum Poisson bracket in the usual way, 

{A, B} := 1 [A, B] = ^ (AB - BA) , 

in in 

then the above basic quantities satisfy the rules identical with the classical ones 

(ESD-dMi), HMD. 

The same concerns the co-moving representants based on differential opera- 
tors H^, ~E A b- The corresponding formally Hermitian operators 

;i ft „ d £i A ft n 9 

v A = -n A = - v a A —, £ s = — e b = ~<p b~q ~ j 

l l OX l l Op A 



and 

tA ■ l = a f L -r ^ L = i\. L - ' L 



J K r = a A 'Pr + S A ' r = A K r, + 



are quantum generators of GAf(iV), GL(U). 

Operators of angular momenta are given by the doubled g-skew-symmetric 
parts of affine momenta: 

S"fa = S a f, — g ac gbd^ d c, L a f) = A a ;, — g ac gbdA d c , 

and 

J a b = J a b — g ac gbd2 d c = L% + S%. 

They are respectively spin, orbital, and the total angular momentum. Similarly, 
the quantum vorticity operator is given by 

\ T A s^A „ac„ 

V b = Zj B — r) T)BD C- 

Canonical linear momentum conjugate to tp 1 a is on the quantum level repre- 
sented by the operator 

pA = h d 

i dip 1 a 

Important: it is not formally Hermitian in L 2 (Q,a), L 2 (Qj nt ,A). Indeed, clas- 
sically it generates additive translations in L(U, V): 

<P 1 A H-> P 1 A +CA- 
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And those do not preserve the Haar measures a, A. But they preserve the 
Lebesgue measures a, I, therefore, P A i is formally Hermitian in L 2 (Q,a) and 
L 2 (Qint, I)- Because of this the Hilbert spaces are more convenient for describing 
quantization of models based on the usual d'Alembert principle. 

However, in all kinds of models the "non- usual" Hilbert spaces may be also 
applied, simply the definition of some operators must be modified. For example, 
affine models may be as well quantized in L 2 (Q,a) and L 2 (Qi nt ,Z), but the 
operators A, B in (|75|) . ([75)) are to be replaced by A', A' given by 

A' := dctLL4p (,l+1)/2 A, B' := dctB-" /2 B. 

More explicitly, 

(A'*)(x jV j) = detL[A]~( n+1 V 2 V{A{x),L[A}<p), 
(B'V)(x,ip) = detB- ,l/2 ^ (x, tpB) . 

Similarly, when considering only the action of A € GL(V) on internal degrees 
of freedom, we define 

A' := dctyr™ /2 A, 

i.e., 

(A'*) (x, ip) = detyr ,l/2 1< (x, Asp) . 

Due to the multiplicative terms, the above operators are unitary in Hilbert 
spaces based on the Lebesgue measures. Their infinitesimal generators are then 
modified by additive correction terms due to which they become formally Her- 
mitian in L 2 (Q, a) and L 2 (Q int , I). For example, T, a b , K a b , S A B are respectively 
replaced by 

'£% = £% + ^n, 'A\ = A\ + ^5\, >V A B = V A B + ^5 A B . 

2% 2i 2i 

It is easy to see that the linear momentum P a , spin S a b, vorticity V B re- 
main unchanged. The finite actions generated by them preserve the Lebesgue 
measures. 

Similarly, when quantizing the d'Alembert models with the use of non-usual 
for them L 2 (Q, a) and L 2 (<5i nt , A), we would have to modify P A i = (h/i)d/d(p l A, 
but we shall not do this here. Nevertheless, it must be stressed that quantiza- 
tions in terms of "non- usual" Hilbert spaces may be convenient when one is 
interested in comparison between various models. 

As mentioned, the best tool when quantizing the d'Alembert model is offered 
by the geometry of Hilbert spaces L 2 (Q, a) and L 2 (Qi nt , I). Then the quantized 
version of thee kinetic energy (|33p is given by the operator 

T = T tr + T lnt - -W'p,p., + I J A BP A t P B 3 g lJ - (77) 
2m 2 

Explicitly, this is a kind of "Laplace operator" in the n(n + l)-dimensional 
Euclidean space: 

h 2 d 2 h 2 ~ d 2 
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Unfortunately, geodetic (potential free) models are non-physical because 
they predict only escaping, non-bounded classical motion and the purely contin- 
uous positive spectrum after quantization (no bounded states). And for realistic 
potentials the variables (x 1 , tp % a) are rather non-adequate. The classical expres- 
sions are non-convenient for quantization because they suffer from 
the embarrassing problem of the ordering of operators. 

There are no such problems with models based on the affine invariance, i.e., 
(p{§|) , (HO]). Let us remind that the first of them is affmely invariant in the 
physical space and isometrics-invariant in the material space. On the contrary, 
the other one is isometries-invariant in the physical space and affinely invariant 
in the body. There are no ordering problems and the quantum operators of the 
kinetic energy may be immediately obtained via the simple replacement of the 
classical linear momentum and affine spin by the operators just written down. 
Therefore, for the quantized versions of (|39f . (|40| we respectively obtain 



= VaVb + -C a cS b S d 

a d2 h2 ^i 1 i d k 9 
2m ox l oxi 2 dtp J a otp 1 B 

where, as previously, C denotes the inverse Cauchy deformation tensor. The 
ordering of operators tp, 8/ dtp is essential, therefore, there appear first-order 
differential operators, respectively, 

^ 2 rB A ,k ® hi rB A 

2 A ° B ~dWc = ~^ A ° B ' 

~~2 A ~dp^A = -^ Ul ^ 1 - 

The second-order terms are obvious: 

h 2 ~ FP- fr 2 ~ fl 2 

rB D ,k , J ° a rr>j I ,J ,k ° 
-Z-l- A C<P BP> D h r, , i 7Tl< - i W AP> B- 



2 0tp K AOtp l c 2 Otp^A<Jtp l B 

There "curved" structure is obvious. Geometrically this is due to the fact that 
the metric tensors on Q given by (|57|) , (|38[) define there essentially Riemannian 
structures with non-vanishing curvature tensors. Unlike this, the d'Alembert 
model is based on the evidently flat, Euclidean geometry with the metric tensor 
(|34|) . All this has to do with strong nonlinearity encoded in geodetic terms 
of classical affine Hamiltonians. And this strong nonlinearity follows from the 
"large" group of assumed symmetries. 

As mentioned in the classical part, there are good reasons to concentrate the 
attention on those metric tensors on Q = Mx Ll(U, V) (those models of kinetic 
energy) which are: 
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1. affinely- invariant in the physical space and simultaneously isometry- inva- 
riant (homogeneous and isotropic) in the material space, 

2. and conversely homogeneous and isotropic in the physical space (isometry- 
invariant) and simultaneously amnely-invariant in the material space. 

It is impossible to satisfy simultaneously both demands [T] and [2] However, if 
translational degrees of freedom are neglected, there exist metrics on Qint = 
LI(£7, V) affinely (or rather centro-affinely) invariant both in the space and in 
the body. They are always pseudo-Riemannian, i.e., have the non-definite hy- 
perbolic signature. Obviously, ([75]) . (|79[) are formally Hermitian in L 2 (Q,a). 
Just like dTTJ) is so in L 2 (Q, a). 

The operators of translational kinetic energy T tr are exactly like in general 
models ([78]) . ([79]) . so we concentrate on the internal parts Ti nt . And they are 
just the very special cases of those in the general formulas (|78|) . (T79|) . Due to 
their very peculiar features it is instructive to quote them explicitly. 

So, for internal degrees of freedom the quantized versions of P9"]) . (|5"0|) are ob- 
tained by the literal substitution of £, £ instead of their classical counterparts, 
so, respectively, 

T int = i^^^EVsV + is^S^ + isVE^, (80) 
21 2A 2B 

T int = -ijft/^E'i + ^E't + -IeSe'i (81) 
21 2A 2B 

with the same meaning of inertial constants /, A, B like previously, ([1]). And 
again the second terms of both expressions are identical; the same is true of the 
third ones. Let us write explicitly 



;£ S A = = -h 2 v k BV l A ——— - inhip k 



d<p k Ad<p l B dip k A 

The geometric interpretation of the affine spin, usual spin, and vorticity as 
generators of transformation groups implies that many quantum expressions 
involving them may be, as it was just seen, directly obtained from classical 
formulas by simple substitution of appropriate operators instead of the corre- 
sponding classical phase-space quantities, so that the difficult ordering problems 
are avoided. 

For example, the very convenient classical expressions (I53[) . (|54[) remain valid 
on the operator level as alternative expressions respectively for (f8TJ|) , (f8"Tj) free 
of the inconvenient transposition term (one with /) : 

Tint^^ + Ac^ + i-IIVlP (82) 

(affine-metrical model), and 

T Jrt = i0(2) + ic(l)» + i I ||S||» (83) 
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(mctrical-amnc model), where a, (3, /i are the same constants as (j52")l from 
classical formulas, and the operator Casimirs C(fc), ||V|| 2 , ||S|| 2 are given by 

C(l) - X A A = £ a Q , C(2) = X A B X B A = £ a fc £ b a , 
and analogously for k > 2 (till k — n), and 

l|V|| 2 = - l -\ A B V B A , ||S|| 2 - -\s\s\. 

The more so the amne-amne model (|55p retains its structure when quantized 
(transformed to the operator form) : 

And similarly, the splitting of the kinetic energy into incompressible and dilata- 
tional parts survives smoothly the quantization procedure. Decomposition (|28|) , 
(|29|) of the affine spin into the spin-shear and dilatation parts has the following 
form: 

V a b = s\ + -p6 a b , V A B = s A B + -p5 A Bl 

n n 

where the trace-less parts are given by 

S a b = ^ a fc ^ d dS a b, S A B = ^ A B Ti D B S A B . 

n n 

They are formally Hermitian operators generating the unitary actions of SL(V), 
SL(U) on L 2 (Q in t, A) in the sense of ([75]). 

A G SL(V) : (A*) (<p) = V{A<p), 
B G SL([/) : (BW) {ip) = V(tpB). 

Incidentally, they are also formally Hermitian in L 2 (Qi nt ,l), because the above 
actions of SL(V), SL(U) are there unitary. 
The operator p has the following form: 



_ H d 
A i dq 



It is formally Hermitian in L 2 (Qi nt , A) (but not in L 2 (Qj nt , I)) and generates the 
one-parameter unitary group of dilatations. This is the group induced by the 
additive translations of logarithmic deformation invariants q a , in particular, by 
the additive translations of their "centre of mass" q. This is the group which 
acts as follows: 

A = e { 6l + : (A*) (ip) = ^(X(f), q^q + £. 

Let us denote the second-order Casimir operator for SL(n, M) by CgL(n)(2): 

CsL(n)(2) = S a bS b a = S A B S B A , 
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and similarly for higher-order ones C SL („)(fc) (but of course CgL(n)(l) = 0; f° r 
orthogonal groups all of the odd orders vanish) . 

On the quantized level the structure of affine-affine, affine-metrical, metrical- 
affine models (respectively, ([ST)) . (j55|) . (|59p) beautifully survives in operator 
language. Namely, one obtains, respectively, 



- 2l C SL(n)(2) + nB) p 2 - T° h + T° u , 

1 1 o I 



Tint " 2(7TA) CsL(n) (2) + 2n(J + A + nB) p2 ^ 2(7* - A*)" * 11 
Tint = 2(7TA) CsL(n) (2) + MI + A + nB) p2+ W 2 A 2 ) M ^ 



Obviously, 



d 2 

" 2 - - h w 



All these formulas are automatically obtained from the corresponding classical 
expressions (|57j) . (|58|) . (|59j) by the formal substitution of operators instead of 
phase-space quantities. It is so because one deals here with generators of the 
underlying transformation groups, quantities of profound geometric interpreta- 
tion. 

Just as in the classical case the quantum unbounded dilatational motion 
should be stabilized by some potential V(q) if the model is to describe quantum 
elastic vibrations. The Hamilton operator splits then into two independent 
mutually commuting parts 



H H s h + H 



clil- 



The same is true for more general doubly isotropic potentials separating explic- 
itly the shape and dilatation dynamics: 



V = V sh (...,q i -<?,...)+V dil (q). 



We have then 



1 

2A~ 

1 



H sh — — C SL ( n )(2) + Vsh,, 



Hsh = 2(7Ta) Csl( " )(2) + 2(P-a*) 11 v|1 +ysh ' 

Hsh = 2(7TA) CsL( " )(2) + 2(P - tf) l|S " 2 + V ° h 

respectively for the affine-affine, affine-metrical, and metrical-affine models. The 
first two of them reduce to the third one when we put 1 = 0. 
And obviously 

1 h 2 d 2 

Hdil = o 71 a i mP 2 + ^dil(g) = ~ Tf i a i m «~2 + v Ml)- 

2n(I + A + nB) 2n{I + A + nB) oq l 
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Van ma y be chosen as some qualitatively satisfactory phenomenological model, 
e.g., Vdii = (k/2)<7 2 , the finite or infinite potential well, etc. 

When V^h = 0, we are dealing with purely geodetic afhnely-invariant Hamil- 
tonians built entirely of the group generators. In such situations one can expect 
solutions of the eigenproblem based completely on some purely algebraic ladder 
procedure. 

Finally, it is interesting to express everything in terms of the two-polar de- 
composition. And now an unpleasant surprise, namely, the automatic replacing 
of classic quantities by seemingly natural operators does not work any longer. 
Namely, although p may be automatically substituted by p = (H/i)d/dq, it is 
not the case with p a , they are not replace by (h/i)d/dq a and J2 a Pa m <P>3"j) is 
not "quantized" to 

-h 2 A{q a ]=-h 2 Y 

1 1 ^ dq a2 

a 

The point is that the additive translations of logarithmic deformation invari- 
ants are not geometrically fundamental operations. So, whereas (|82|). (|83|) . (I84|) 
are automatically (— fi 2 /2)-multiples of the corresponding Laplace-Beltrami op- 
erators, there is no such automatism with the two-polar expression of these 
operators. Fortunately, there are no problems with the spin and vorticity op- 
erators S 4 j, ~V A b and with operators r a t, t a f, corresponding to the classical 
quantities (|32p . The reason is again their group-theoretical interpretation: spin 
and vorticity generate respectively spatial and material rotations. And the op- 
erators r%, t% are their representations in terms of the principal axes of the 
Cauchy and Green deformation tensors, 

r^b = L a iL\S l j, t a b — — R a aR B bV A b\ 

the ordering of operators just as written here. Just as in the classical theory, r a b, 
t a b are generators (in the quantum-Poisson-bracket sense) of the right action of 
SO(n, E) on the quantities L : 1" -> V, R : 1" -> U 

L^LU, R^RU, UeSO{n,M). 

Just as r a b, t a b, the operators S 4 j, V A b act only on generalized coordinates x^, 
parameterizing respectively L and R (some Euler angles, rotation vectors, 
first-kind canonical coordinates, and so on). Any of the mentioned classical 
quantities S l j, V A b, r a &, t% has the following form: 

r(x) P (x)», $"(i/)p(y) M) 

where p{x) lll p{y)fi are canonical momenta conjugate respectively to x^, y^. Due 
to the group-theoretical meaning of the mentioned quantities, the corresponding 
quantum operators are given respectively by the operators 

V(x)— -aHv) — ■ 
i 1 1 W' i 9 (V) dy»' 

the ordering just as explicitly written. 
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In analogy to (|62[) we introduce the operators 



M = -r-t, N = r-t. 
The kinetic energy operator for affine-affine model is then given as 

h 2 _ h 2 B d 2 



rpiaff — aff ■ " t~v 

1 int - - o~7 



2A 2A{A + nB)dq 2 

1 v M 2 ab 1 v N 2 ab 

o,6 b 2 o,6 2 



where 



A ~p7^d^ A 9^^^5^2 + ^ dq a !hf ( ' 

and P A is given by l[7IJ). The expression ([55]) exactly equals -(?i 2 /2)A(r ), 
where the Laplace-Beltrami operator A(r°) is built of the configuration metric 
r° (j4T?|) , i.e., corresponds to the classical expression (|4"5|) , (|45[) . It is seen that T>\ 
differs from the R n -Laplace operator ^ Q d 2 /dq a2 by some first-order differential 
operator. This is just the mentioned breakdown of the naive classical analogy 
between p a and (h/i)d/dq a . The reason of this breakdown is that the additive 
translations 

q i— * q + u 

do not preserve the measures A, a. Because of this their argument-wise action 
on wave functions is not unitary in L 2 (Q, a), L 2 (Qi n t, A). Incidentally, it is not 
unitary in L 2 (Q,a), L 2 (Qint,0 either. And infinitesimal generators (h/i)d/dq a , 
(H/i)d/dQ a are not formally self-adjoint in these Hilbert spaces. 

The affine-metric and metric-affine models are respectively given by 



rriaff— met i\ 
T int = - o~~ 



ft 2 

2a" X ~2f3dq 2 



_^ T ^ 1 y- N « b | 1 || V |. a 

32 «^Bh 2 ^ 32a f-ch 2 ^ 2/J 1 ' 



fc2 fc2 a2 

rpmet-aff _ _ _T) _ 



2a 2/3 dq 2 

1 v M 2 b 1 ^ Ng b 1 

32a f^sh 2 ^ 32a ^ch 2 ^ 1 2 ^ 



Quite similarly, quantizing the doubly-isotropic d'Alembert model we obtain 



T5L A = --D, + 



im 2/ ' 8/ (Q a - Q b ) 2 81 {Q a + Q b ) 2 : 
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where 



d 



Pi- 



d 



Pi ^ dQ a dQ a 



- V— — V 

~ ^ dQ a2 + ^ 



dlnPi d 
dQ a dCF*' 



(87) 



and the weight factor Pi is given by (|72[). 

The ordering of non-commuting operators Da, D; is just as explicitly written 
here. There are no other ordering problems because the operators M a fc, N a £, 
(equivalently f a h, t a b) do commute with deformation invariants q a , Q a . 

Let us observe that the first-order differential operators in (|86|) . (|87f may be 
eliminated by introducing modified amplitudes ip given respectively by 



(88) 



Then the action of D-operators on if? is represented by the action of operators 
D given respectively by 



2A 



D 



A = 



2 A 



E 



~,a2 



21 



; D, 



21 



E 



d 2 



dQ 



it 2 



Vi, 



where V\, Vi are auxiliary "artificial" potentials 
Vx = 



n_\_ h 2 l x ^ ( dP A 

2AP 2+ lA~Px^ [ihf) 
tf-_l_^fdPi 



X 

h l 
'Yip? 



\dQ a 



In other words 



Dip = \/PD*. 



All other terms of the kinetic energy operators commute with v P interpreted as 
a position-type operator. Obviously, the same concerns usual potential terms V. 
Therefore, finally, the Hamilton operator H = T + V is represented in ip-terms 

b y H . 

Hip = VPH*, 

where analytically the action of H differs from the action of H in that the 
D-operators are replaced by the usual IR™-Laplace operators and additional 
V-potentials appear and are combined with the "true" potentials V. The sta- 
tionary Schrodinger equation, i.e., eigenequation 



is equivalent to 



= E^f? 



Hip = Eip. 



This eigenproblem is meant in the Hilbert space based on the modified scalar 
product without the weight factor P in the integration element, 

(fil^) = / Tpl^dq 1 ■ ■ ■ dq n dfj,(L)dfi(R) 
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in affine models, and 

(<^iM = J Tpl^dQ 1 ■ ■ ■ dQ n d^{L)d^(R) 
in d'Alembert models. Obviously, 

(piM = (*l|*2>- 

There is no real simplification in replacing 4" by <p because instead of the com- 
plicated first-order differential operator the equally so complicated potential V 
appears. 

In geodetic problems and in problems with the doubly isotropic potentials 
V(q 1 , . . . ,q n ), in particular, with the stabilizing dilatation potentials Vdn(q), the 
above Schrodinger equations may be reduced to ones involving only coordinates 
q 1 ,. . ., q n , because the action of H on the (L, i?)-dependence of wave functions 
may be algebraized. This is based on the generalized Fourier analysis on the 
compact group SO(n,R). 

To simplify the treatment we identify analytically Qmt with GL + (n,R) and 
use the matrix form of the two-polar decomposition <p = LDR -1 . According 
to the Peter- Weyl theorem, the wave functions may be expanded in (L, R)- 
variables with respect to matrix elements of irreducible unitary representations 
of SO(n, M). Their expansion coefficients are functions of deformation invariants 
q a , or equivalently, Q a . Let fi denote the set of irreducible unitary represen- 
tations of SO(n,R) (more precisely, the set of their equivalence classes). Obvi- 
ously, due to the compactness of the group SO(n,R) these representations are 
finite-dimensional; their dimensions will be denoted by N(a). In the physical 
three-dimensional case fl is the set of all non- negative integers s — 0,1,2,... 
and N(s) = 2s + 1. If for som e reasons we replace the rotation groups by their 
universal coverings SO(n, R) and so admit half-integer angular momenta, then 
fi is the set of all non-negative half-integers and integers s = 0, 1/2, 1, 3/2, . . . 
and again N(s) = 2s + 1. Obviously, in the planar case = Z is the set of all 
integers and N(m) = 1 for any m e Z (Abelian group). 

Let D a be N(a) x N(a) matrices of irreducible representations. Then the 
mentioned expansion has the following form: 

N(a) N{f3) 

y(<p) = *(L,D,R)= ]T ]T Y j V a mn {L)f:J i {D)V^ kl {R- 1 ). (89) 

a,/3€nm,n— 1 k,l— 1 

The non-uniqueness of the polar decomposition implies that the deformation 
invariants q 1 , . . . , q n (Q 1 , . . . , Q n ) are very complicated indistinguishable paras- 
tatistical "particles" in R. There is no place here to get into more details. The 
point is that the reduced amplitudes jZk as functions of q 1 , . . . , q n must satisfy 

ml 

certain conditions due to which the resulting $ as a function of L, D, and R 
does not distinguish triplets (L, D, R) representing the same configuration <p, 
i.e., *(Li,£>i,i?i) = y{L 2 ,D 2 ,R 2 ) if LxD^ 1 = L 2 D 2 R^ 1 . This is simply 
the condition for ^ to be a one- valued function on the configuration space Qint- 
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One can consider the matrix elements T> a m k as explicitly known. And in 
fact, they are deeply investigated special functions on the orthogonal groups 
SO(n, K). In the physical case n = 3 they are well-known functions T>^ m k 
found by Wigner. Here j = 0,1,2, ... or, if we replace SO(3, R) by its universal 
covering SU(2) (for the general n, SO(n, M) is replaced by the group Spin(n)), 
j = 0, 1/2, 1, 3/2, . . .. And according to the standard convention m = —j, —j + 
1, . . .,j—l,j, k = —j, + — for the fixed j, m, and k have (2j + 1)- 

element integer range with jumps by one both for the integer and half-integer 

i- 

The operators S z j, V A b, ? a fc, t a f, when acting on functions V a m k may be 
replaced by some standard algebraic operations. This enables one to reduce the 
Schrodinger equation for the wave functions \& depending on n 2 variables tp l A 
to some eigenproblems for the multi-component amplitudes f a/3 depending only 
on the n deformation invariants q a . Therefore, in a sense, the problem may be 
reduced to the Cartan subgroup of diagonal matrices ip (the maximal Abelian 
subgroup in GL(n,R)). 

In geodetic models and in models with doubly isotropic potentials (ones de- 
pending only on deformation invariants; dilatations-stabilizing potentials V(q) 
provide the simplest example), the labels m, n in (|89p are good quantum num- 
bers. The Hamilton operator H commutes with the operators of spin and vor- 
ticity, i.e., S z j, V A B . Also the representation labels a, (3 £ f2 are good quantum 
numbers. They are equivalent to the systems of eigenvalues of the Casimir 
invariants built of S, V: 

C(S,p) = S l k S k m ---S r z S z t , C{V, P )=V A kV k m---V r zV z a (90) 

(p factors). These eigenvalues will be denoted respectively by C a (p), C^(p). 
Obviously 

C(S,p) = C(r,p), C(V,p) = c(t,p). 

The above Casimir invariants vanish trivially for the odd values of p; so in the 
physical case n = 3 there is only one possibility: C(S, 2), C (V, 2). Due to the 
peculiarity of dimension three, where skew-symmetric tensors may be identified 
with axial vectors, it is more convenient to use 

||S|| 2 = -^S\S b a , ||V|| 2 = -lv A B v B Al 

i.e., (-l/2)-multiples of C(S, 2), C (V, 2). The point is that for n = 3 
||S|| 2 = Si + S 2 + S 2 , ||V|| 2 = V 2 + V 2 + V 2 , 

where 

S a = 2 £ afc c S b c , Va = 2 £a b C V B c- 

The raising and lowering of indices is meant here in the sense of orthonormal 
coordinates (Kronecker-delta trivial operation). The same convention is used 
for f° 6 , t a fc , i.e., 

f = -f h c r b t = -£ h c t b 

1 a n c ao L c: L a — n c ao L c- 
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Obviously, 

||?|| 2 = ||Sf, ||tf = ||V|| 2 . 
The corresponding eigenvalues are given by 

C(s, 2) = h 2 s(s + 1), C(j, 2) = H 2 j(j + 1), 

where s,j are non- negative integers or non-negative integers and positive half- 
integers when GL + (3,R), SL(3,R) are replaced by their coverings GL + (3,M), 
SU(2). 

It is convenient to use multi-component wave functions with values in the 
space of complex N(a) x N({3) matrices ((2s + 1) x (2j + 1) in the physical case 
n = 3): 

*(y>) = * a/3 (L, D, R) = V a (L) / Q/3 (D)V^ (R) , (91) 

where f a/3 are complex N(a) x N((3) matrices-reduced wave amplitudes de- 
pending only on the deformation invariants. In the physical three-dimensional 
case, when T> a mn are Wigner special functions T> s mn , we, as usual, take m,n 
running from — s to s and jumping by one (also in the "spinorial" case when s 
is half-integer). And then 

||S|| 2 * SJ = h 2 s{s + i)$ sj , ||v|| 2 ^' = h 2 j(j + i)t> sj , 

And similarly, when the values n, k in the superposition (|89|) are kept fixed and 
wc retain only the corresponding single term, for the resulting ^ we have 

Let us now describe in a few words the aforc-mcntioncd algcbraization proce- 
dure in the sector of (L, i?)-degrees of freedom. If a compact group counterpart 
of the usual Fourier-transform algebraization, where d/dx a is represented by 
the point-wise multiplication of the Fourier transform of f(x), f (fc) by ik a . 

Let us introduce some auxiliary symbols. 

The group SO(n, R) may be parameterized by the first-kind canonical coor- 
dinates w, namely, 

W(u) = exp Qo^O , 
where the basic matrices of the Lie algebra SO(n, M)' are given by 

(E b a ) d = S b dS C a ~ S bC 5 a d, 

and the matrix u> is skew-symmetric in the " cosmetic" Kronecker sense. There- 
fore, independent coordinates may be chosen as uj a bi a < b, or conversely. How- 
ever, for the symmetry reasons it is more convenient to use the representation 
with the summation extended over all possible uj a b- 
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To be more "sophisticated" , the groups SO(V, g), SO(U, r)) are parameterized 
as follows: 

W(cj) = exp Qw^'i) , W(u>) = cxp (^Q A B E B J 
where E l j, E a b are basic matrices of Lie algebras SO(V, g)', SO(U, rj)' given by 

(E i j )\ = Wi-g^gn, (e A b) C d = S A D 5 C B - r, Ac r, BD . 

The skew-symmetry of u) in the above exponential formulas is meant respectively 
as follows: 

LU a b = -g ac g bd uj d c , Cj a b = -7] AC r]BD^ D C- 
Now matrices of irreducible representations T> a are given by 

V a (L(0) - cxp (^r b M ab )j , V a (R(r)) = cxp (^r\M ab }j , 

where I and r denote the w-parameters, respectively, for the L- and i?-factors of 
the two-polar decomposition. The anti-Hermitian matrices M a will be expressed 
by the Hermitian ones S a , 

S aa b = -M aa b . 
1 

The commutation rules for M aa b are expressed through the structure constants 
of SO(n,R), 

[M s abl M s cd ] = -g ad M s cb + g cb M s ad - g bd M s ac + g ac M s bd , 
and therefore 

— [S j ab, S j cd] = g a dS j cb - g C bS j ad + 9bdS 3 ac - 9acS j b d- 

in 

Indices here are shifted with the help of g ab ; as a rule we use orthonormal 
coordinates when g ab = 5 ab . 

In the physical three-dimensional case when we put 

qj lp bcqj 

^ a ■ 2 ° bc^ 

we obviously have 

— \S 3 a , S J '(,] = e ab c S J c . 
in 

From the fact that V a are representations and (i/h)S k i, (i/h)V A B, (i/h)r a b , 
(i/H)t a b are infinitesimal generators of left and right orthogonal actions on the 
(L, i?)-variables it follows immediately that 

= gd^ap^ r a b y af3 = V a {L)S aa b f af3 {D)V r3 (iT 1 ) , 

V A B *f al3 = ■$ af3 S f3A B , t a b ^ a(3 = V a {L)f af3 (D)S t3a b V f3 (R- 1 ) . 



■59 



Therefore, spin and vorticity act on the wave amplitudes ^> a/3 as a whole, and 
in a purely algebraic way. On the other hand, to describe in an algebraic 
way the action of r a b, t a b, one must extract from 1 i ,af3 the reduced amplitudes 
f a ^(q l , . . . , q n ). And it is only this amplitude that is affected by the action of 
r a b, t a f, according to the following rules: 



f*fi S aa b f afl , t a b : f af3 f^S fia b. 
It is very convenient to use the following notation: 

gda ^j-afl S aa bf a ^ a bf a ^ f a ^ 

As T> a are irreducible, the matrices 

/~i(qcl ~.\ . Qaa nab Qau Qaw 

^ W ? Pi ■ — a b^> c ' ' ' <J a 

(p factors) are proportional to the N(a) x N(a) identity matrix, 

C(S a , P ) = C a (p)I N(a) , 



where C a {p) are eigenvalues of (|90| . 

In particular, in the physical case n = 3 we have 

||S|| 2 *^' = \\r\\ 2 ^ = h 2 s(s + 1)W, 

||V|| 2 * SJ ' = ||t|| 2 ^' = h 2 j{j + i)y sj , 



S„>1'-' .s-...>l'-'. 

V ■ ' = V s & a, 



where S s a are standard Wigner matrices of the angular momentum with the 
squared magnitude h 2 s(s + 1). Multiplying them by (i/h) we obtain standard 
bases of irreducible representations of the Lie algebra SO(3,M)'. For the stan- 
dard Wigner representation the following is also true: 



ml i 



v 3 *Z = raw 



ml ' 



Similarly, the action of r, t operators is represented by the following operations 
on the reduced amplitudes: 



f sj _ f s 3 gJ a = gJJS 



In particular, 



r 3 



J ml 



hm fZ 



•> ml 



S3 



Using the well-known orthogonality relations for the matrix elements of irre- 
ducible unitary representations V a mn [3TJE5] we can rewrite the scalar product 
in the following form: 



i 



a, pen 



N(a)N(P) 



Tr 



(ff + {<f)ff{<l b ))Pxdq l ---dq n , (92) 
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where Pa is the weight factor given by (fTTj) . and the argument symbols like q a 
are abbreviations for the system (q , . . . , q n ). The trace operation is meant in 
the sense of matrix two-indices: 

N(a) N(fi) 

Tr (ff + ff) = E E A t\ h tv 

n,m—l k,l—l 

If no superposition over m, I in (|89| is performed and we use the matrix- valued 
wave functions (|9ip . the trace operation is meant in the usual sense. 

Obviously, if we use the modified wave functions tp (|88p . then the scalar 
product expression is free of the weight factor P\, 

<*>. w - E n ^^m / ^ (< + w> 9 f fcfl) 

where, obviously, g = \fP\f '. Quite analogous formulas are true for the d'Alem- 
bert models; simply P\ is replaced then by p ([72]) and instead g a we use their 
exponential functions Q a , 

<*i|*»> = E o jV(a )W) / ^ ^ (Qb) ) PldQl ' ' ' 

= E^ jy^cg) / Tr ( g "" + (Qa) dQl ■ ■ ■ dQn - 

Remark: it is implicit assumed in the above formulas that the Haar measure 
on the (L, P)-manifolds is normalized to unity (the total "volume" of the cor- 
responding manifolds equals one). 

We restrict ourselves to Hamiltonians of the form H = T + V with some 
doubly-isotropic potentials V(q 1 , . . . , q n ), in particular, with some dilatation- 
stabilizing potentials V(q) (affinely-invariant geodetic incompressible models). 
The energy eigenproblem, i.e., stationary Schrodinger equation 

H* = P* 

is equivalent to the infinite sequence of eigenequations for the reduced multi- 
component amplitudes f a @: 

The simultaneous spatial and material isotropy imply in the N(a) x iV(/3)-fold 
degeneracy, i.e., for every component of the N(a) x iV(/3)-matrix amplitude f a " 
there exists an N(a) x iV(/3)-dimensional subspace of solutions, just as seen from 
the symbol used in 

ml 

The reduced Hamiltonians 
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are TV (a) x TV(/3) matrices built of second-order differential operators. 
For the affinc-affine model of the kinetic energy we have 



f, 2 1 (s Pa b - P a b ) 

T af3 fa f3 = _ 'J_r>f Q ' 3 4. _L_ V ^ '—f 

J OA J 39 4 „- 2q°--q* 



2 

a/3 



2A 32Aj^ sh 2 

1 (^ + ^V) 2 ^ ,2 

32 A ^ ch 2 ^^ 2A(A + nB) dq 2 3 

For the spatially metrical and materially affine model we obtain 



, 2 



T a/3 / a/3 = _ ^_ D/ a/3 + J^, a2 )/^ + — V 
J On 1 On V ' > J 10rv L^i 



fa/3 



2a 2^ v ' /J 32a ^ sh 2 ^ 



_ _V-^ J_ faf 3 _ ^_^_ faf 3 m) 

32a ^ ch 2 2l^ 7 2/3 9g 2 7 ' 1 j 

0,6 Ln 2 M y 

where C(a, 2) is the a-th eigenvalue of the rotational Casimir ||S|| 2 , thus, 

~S ai j S a ' i = C{a,2)l N(a) . 

Obviously for the physical dimension n — 3, f a/3 — f^, we have C(s,2) = 
h 2 s(s + 1). And similarly for the spatially affine and materially metrical model 
we have 



h 2 1 1 (^b-^b) 2 



•a/3 



1 ^ f ^% + S*<0 ft 2 £2 

32a ^ C h 2 2^1 7 2/3 9 9 2 7 ' 1 j 

a, 6 2 

where C(/3, 2) appears as the /3-th eigenvalue of the vorticity Casimir ||V|| 2 , 
and just as previously for n = 3, / Q/3 = / s -?, we have C(j,2) = h 2 j(j + 1). 
It is so as if the doubly affine background (T affinely-invariant in the physical 
and material space) was responsible for some fundamental part of the spectra, 
perturbated by some internal rotations of the body itself or of the deformation 
axes. This perturbation and the resulting splitting of energy levels becomes 
remarkable when (i is small, i.e., when the inertial constants /, A differ slightly. 
The suggestive terms 

2^ S(S + 1) ' 2^ + 1) 
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as contributions to energy levels are very interesting and seem to be supported 
by experimental data in various ranges of physical phenomena. 

Finally, let us quote the corresponding form of T Q/3 for the quantized d' Alem- 
bert model: 

T a/3 fCt /3 = _ fa , J_ V }_ fa/3 

J 21 11 ^81^ (Q a ~Q b ) 2 J 

a,b 



1 + 
8/ 2-? (Q a + Q h ) 2 1 



In this way the problem has been successfully reduced from n 2 internal degrees 
of freedom (physically 9, sometimes 4) to the n purely deformative degrees of 
freedom (physically 3, sometimes 2). The price one pays for that is the use 
of multi-component wave functions subject to the strange parastatistical condi- 
tions in the reduced q a -variables. The particular values of labels a, (3 and the 
corresponding matrices S a a b, a b describe the influence of quantized rotational 
degrees of freedom on the quantized dynamics of deformation invariants. It is 
interesting that on the classical level there is no simple way to perform such a 
dynamical reduction to deformation invariants. 

For any reduced problem with a, (3 labels the quantity 

pi?) :=Tr(ff + {q«).ff (q h ))P X (q c ) 

is the probability density for finding the object in the state of deformation 
invariants (q 1 , . . . , q n ). More precisely, p(g 1 , . . . , q^dq 1 ■ ■ ■ dq n is the probability 
that the values of deformation invariants will be detected in the infinitesimal 
range dq 1 ■ ■ ■ dq n about the values (q 1 , . . . , q n ). 
Similarly, performing the integration 

p(L,R) = J * (L; q a ;R)y (L; q a ; R) P x (q) dq 1 ■ ■ ■ dq n 

one obtains the probability density for detecting the "gyroscopic" degrees of 
freedom L, R (equivalently, the Cauchy and Green deformation tensors) in some 
range of the configuration space. Obviously, this distribution is meant in the 
sense of the Haar measure fi. The integrals 

V a l = fit? (q a ) f t (q b ) Px (ql dq 1 ■ ■ ■ dq n 

ml J ml ml 

are probabilities of detecting the particular indicated values of angular mo- 
menta and vorticities, C(a, 2), C(/3, 2), hm, hi, hn, hk. In the physical three- 
dimensional case they are respectively h 2 s(s + 1), h 2 j(j + 1), hm, hi, hn, hk, 
where s, j are non-negative integers, m,n — — s, . . . , s, k,l = —j, . . . , j, jumping 
by one. Particularly interesting are 

Pml — P ^ 

n,k 
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because they refer to the constants of motion S*j, V A b and to "good" quan- 
tum numbers a,(3,m,l. Except the special case n = 2, r a b, t a & are not con- 
stants of motion, and k,l are not "good" quantum numbers. Quite analogous 
statements are true for the quantized d'Alembert model; the only formal dif- 
ference is that the integration element in the manifold of invariants is given by 
P l {Q\...,Q n )dQ 1 ---dQ n . 

Let us now write down the explicit formulas for the physical three-dimen- 
sional case. For the affme-affine model (1951) we have now 



fc2 t2 r> £)2 

aS ~ aSj 2A XJ 2A(A + 3B)dq 2j 



1 3 

— T 



a=l blL 2 

1 ^ (S\) 2 f s * + 2S\f°*& a + f a *(& a ) 2 
16 A ^ C h 2 ^ ■ [M) 

where S J a are the standard Wigner matrices for j-angular momentum, i.e., 
h 2 j(j + l)-magnitude, and for any a-th term of both summation we have obvi- 
ously b 7^ a, c 7^ a, b ^ c. Obviously, it does not matter in what an ordering q b , 
q c are written, because the denominators are sign-non-sensitive. Da is given by 
((86)) , where explicitly 



P A = |sh (q 2 - q 3 ) sh (q 3 - q 1 ) sh (q 1 - q 2 ) \ . 

And similarly, using the abbreviated form, we can write for the metrical- affine 
(|94|) and amne-metrical ((95)) models, respectively, as follows: 



TSet-aff = T^_ aff [A » I + A] + 2(/2 _ A2) h 2 s{s + 1), 

T&_met = T^_ aff [A » I + A] + 2(/2 ^ A2) fi 2 + 1), 

where, obviously, T a g_ aff [A h> / + A] is obtained from Tfi ~ simply by 
replacing A with a = I + A. 

The doubly-isotropic d'Alembert model in three dimensions has the following 
form: 

/<'„ , 1 ^ [S'affi -2S' a f'lS! u + f'*{Si u ) 2 



21 lJ 4J^J (<3 b - <3 C ) 2 



i_ \- (s a a yf°3 + 2s° a f°>s' a + r(S' a y 

41^ (Q b + Q c ) 2 



a=l 



with the same convention as previously, D; given by (|87p . and explicitly 



{(Q 2 ) 2 - (Q 3 ) 2 ) {(Q 3 ) 2 - (Q 1 ) 2 ) ((Q 1 ) 2 - (Q 2 
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Let us mention that in principle half-integer angular momentum of extended 
objects may be formally introduced by replacing the group GL(3, R) by its uni- 
versal covering GL(3, M). There are some indications that the physical usefulness 
of such models is not excluded. Formally, the procedure is as follows. In (189|) 
specialized to n — 3 we replace the group SO(3,R) by its covering SU(2) and 
write the following expression involving the known Wigner matrices V s mn : 

*(u,Z>,«)=£ E E V s mn {u)fd{D)W k l(v- 1 ), (98) 

s,j m,n— — s k,l — —j 

where u,v e SU(2), D S Diag(3,R) C GL(3,R), and both the integer and half- 
integer values of s, j are admissible. However, if the function of triples (u, D, v) 

is to represent a function on GL + (3,R), then the values of s,j in the above 
series must have the same "halfness", i.e., either both s,j in (|98[) are integers 
or both are non-integers. And no superposition between elements of these two 
function spaces is admitted (a kind of superselection rule). The point is that 
for such " halfness-mixing" superpositions the squared modulus ^"5 would be 
two-valued from the point of view of SO(3,R). This would be violation of the 
probabilistic interpretation of '5 in GL + (3,R). If there is no mixing, then in the 
case of superposing over half- integer s,j in ([98]) the resulting is two- valued 



on GL+(3,R), i.e., it does not project from GL+(3,R) to GL+(3,R) but 
does project, i.e., it is single-valued in GL + (3,R). 

The simplest possible situation in ([96]) . ([97]) is s = j = 0, i.e., purely scalar 
amplitude / 00 . Then T 00 reduces respectively to 

nn ft 2 ^ d 2 m ft* 

2A X 2A{A + 3B)dq 2 ' ~ 21 U 

i.e., there is no direct contribution from internal degrees of freedom. 

If we admit half-integers, then the next simple situation is s = j = 1/2. Then 

SVa o = (h/2)a a , where a a are Pauli matrices. Therefore, (S^ 2 a ) = {h 2 /4)l 2 . 

Finally, let us briefly describe the two-dimensional situation, i.e., "Flat- 
land" p], n = 2. Obviously, it may have some direct physical applications 
when we deal with flat molecules or other structural elements. But besides, 
the two-dimensional models shed some light on the general situation and en- 
able one to make it more comprehensible and lucid. Indeed, let us observe 
that the expressions ([65]) and ([68]) (without the last p 2 -term) are superpositions 
of two-dimensional clusters corresponding to all possible R 2 -subspaces in R™. 
Obviously, these terms in general are non-disjoint and for n = 3 they simply 
cannot be disjoint (all two-dimensional linear subspaces in R 3 have intersec- 
tions of dimension higher than null; if different, they always intersect along 
one-dimensional linear subspaces). 

There are some very exceptional features of the dimension n — 2. They are 
very peculiar, in a sense pathological. But nevertheless the resulting simplifica- 
tions generate some ideas and hypotheses concerning the general dimension. Of 
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course, later on they must be verified on the independent basis. Let us begin 
with the classical description. 

The one-dimensional group of planar rotations SO(2, R) is Abelian, therefore, 
p = p = S, t = t = — V. In doubly-isotropic models S and V are constants 
of motion and so are jo, r, M, N if n = 2. It is not the case for n > 2, where, 
as always in isotropic models, S, V are constants of motion but p, f do not 
equal S, —V and are non-constant. But it is exactly the use of p, r and their 
combinations M , N that simplifies the problem and leads to a partial separation 
of variables. In two-dimensional space these things coincide and the problem 
may be effectively reduced to the dynamics of two-deformation invariants both 
on the classical and quantum level. The two-polar decomposition ip = LDR^ 1 
will be parameterized in a standard way; using the matrix language we have 



L = 
D = 



cos a — sin a 
sin a cos a 

Q 1 
Q 2 



R = 



cos (3 — sin /3 
sin (3 cos (3 



cxp q 






cxp q 2 



To separate the dilatational and incompressible motion we introduce new vari- 
ables: 

l/i Q\ 2 1 

1 = o [l + 1 ) ' x = q -q . 



Their conjugate momenta are given by 

p=P!+p 2 , Px 



1 



(P2-Pl) ■ 



Angular velocities are given by the following matrices: 



X 



dL l T ~\dL 
— — L = _L — — 

dt dt 



# = ^ R -i = R -i^ 
dt dt 



da 


' 


-1 


~dt 


1 





d(3 


" 


-1 


~dt 


1 






Spin and vorticity essentially coincide with canonical conjugate momenta p a , 
P/3, i-e., 



S = p = p = p 



1 

-1 



V 



-r = -t = pp 



1 

-1 



With this convention the pairing between velocities and momenta has the form: 



da l m ,„ . 



The diagonalizing quantities M = —p — f , N = p — f are also expressed by 
matrices 



M 



1 

-1 



N = n 



1 

-1 
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where 

m=pp-p a , n=pp+p a . 
In some formulas it is convenient to use modified variables 

7 =I(/3-a), S=\{fi + a). 

Their conjugate momenta just coincide with the above m, n, i.e., 

p 7 = m, p s = n. 

The magnitudes of S, V have the form: 

||S|| = b Q | = l|n-m|, ||y|| = | P/3 | = i|n + m|. 

As mentioned, p a , pp, m, n are constants of motion because in doubly isotropic 
models a, (3 are cyclic variables. The corresponding amnc-affme, metrical-affine, 
and afhne-metrical kinetic energies of internal degrees of freedom are respec- 
tively given by 

T aff-aff _ J_ / 2 , 2\ B 2 

T " 2A^ 1+P2j 2A(A + 2B) P 



+ 



1 m 2 1 



16A sh 2 2f^I 16A ch 2£_ai' 



T met-aff = T aff-aff [A ^ I + A] + _ (n - m) 2 , 

T aff-met = T aff-aff ^ ^ j + ^ + * (n + m) 2 , 

Q{1 A ) 

where, as usual, 7~ aff ~ aff [A I + A] denotes T aff ~ aff with A replaced by I + A. 

Separating dilatational and incompressible motion we obtain respectively 
the following expressions: 

aff-aff P 2 , Vl . (Pa-Pp) 2 (Pa+Pp) 2 



/T-arr— an i , r x , 



inl 4(A + 2B) A 16Ash 2 f L6Ach 2 § 



^mct-aff = T aff-aff / + A] + 



J 2 - A 2 ' 

aft-met /7-aff— aff | 



Jo 2 

^j-att-met = ^att-att J + ^] ■ » 



J 2 - A 2 ' 

Canonical momenta p a , pp, or equivalcntly, m, n, are constants of motion 
and their Poisson brackets with the variables q, x, p, p x . Therefore, if we 
are interested only in the evolution of variables q, x but not in that of a, 
P, we can simply replace p a , pp, m, n in the above expressions by constants 
characterizing a given family of solutions. The are effective coupling constants 
for the interaction between deformation invariants q 1 , q 2 . The sh _2 (a;/2)-term 
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controlled by m is always repulsive and singular at the coincidence x = (non- 
deformed shape but dilatation admitted). The ch _2 (a;/2)-term controlled by n 
is attractive and finite at x = 0. At large "distances" of deformation invariants, 
\x\ — > oo, attraction prevails if and only if |m| > |n|, i.e., if p a , pp have the 
same signs, p a pp > 0. If |m| < |n|, i.e., p a pp < 0, then the time evolution 
of x is unbounded. This is just the very special (n = 2) example of that was 
said formerly, namely that in the incompressible and affinely-invariant geodetic 
regime there exists an open family of bounded motions ("elastic vibrations") 
and an open family of unbounded motions ("dissociation", decay). If the total 
deformative motion is to be bounded, then some dilatations-stabilizing potential 
V(q) must be included into Hamiltonian. But even if there is no x-dependent 
potential, our affine geodetic model in the non-compact configuration space of 
incompressible motion may encode bounded elastic vibrations. The same is true 
for n > 2, however the situation is more complicated then because M a b, N a b are 
not constants of motion and also undergo some vibrations. 

Analogous statements are true on the quantum level. The Haar measure in 
our coordinates is given by 

dX (a; q, x; (5) = \shx\ dqdxdad[3 , 

its weight factor equals Pa = |shx|. The wave functions iff are expanded in the 
double Fourier series: 

ima inf3 



* (a; «,*;/?)= ]T f mn (q, x)e ma e 



This is obviously the Peter- Weyl theorem specialized to the two-dimensional 
torus group T 2 . Our integers m,«eZ are just the labels a, [3 from the general 
theory. 

For the affinc-affinc model the reduced operator of the kinetic energy is given 

by 

fc2 fc2 Q2 fmn 

rymn rmn __T) f mn J 

aB ~ aSj " ' A XJ A(A + 2B) dq 2 
h 2 (n-mf _ h 2 (n + m) 2 
16Ash 2 f 16Ach 2 | 

where now D> is expressed as follows: 

D - — !— — (\shx\ — 
\shx\ dx \ dx 

Similarly for the metric-affine and affine-metric models we obtain respectively 

/ft 2 to 2 

J -met-aff ~ J -aff-aff \ A ^ 1 + A \ + j2 _ ^2 ' 

7ft 2 n 2 

TSf"_ mot = T^"_ aff [A / + A] + j 2 _ ^42 ' 
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To avoid the purely continuous spectrum one must include into Hamiltonian 
at least some dilatations-stabilizing potential V{q). The problems is then (as 
usual) separable in (x, g)-variables. And just as on the classical level, for affinely- 
invariant incompressible dynamics, i.e., for the x-sector of the above operators, 
there exists discrete spectrum if \n + m\ > \n — m|, i.e., if mn > 0. This is the 
quantum bounded motion. 

In three-dimensional problems the above condition will be replaced by some 
more complicated one between quantum numbers labelling the reduced ampli- 
tudes /. 

Let us observe that in more general, not necessarily geodetic, problems in 
two dimensions with explicitly separable potentials V(q,x) = Vdii (<z) + Vshix), 
the Schrodinger equation 

j_jmn jran jmn 

where H m " = T mn + V d n(q) + V sh (x) with any of the above T mn , splits into two 
one-dimensional Schrodinger equations. The reduced wave function is sought in 
the form 

f mn (q,x)^ i p mn (x) X (q), 
where tp mn , \ satisfy the following eigenequations: 

h 2 d 2 y 

Hmn mn ,„mn n ,, A . jr rp 

sh ip = E sh ip , HdilX = _ 4(AT2B) rf? 

Here the shear-rotational Hamiltion operator H^™ is given by 

™ ^ n | h 2 (n-m) 2 h 2 (n + m) 2 ^ 

H sh _ aff _ aff - - 7 Da + 16Ash 2 f 16j4ch 2 f + V *W 

in the affine-affine model, and by 



Ih 2 m 2 „ m „ „ m „ Ih 2 n 2 



Tjmn -rrmn , -wrran rrmn , 

"sh-met-aff — n sh-aff-aff "+" p _ ^2 ' "sh-aff-met — n sh-aff-aff ' j2 _ ^2 

respectively in the mctric-afhne and afhnc-mctric models. Obviously, the total 
energy equals E = E s h + E^n- It is seen that the main point of the analysis is 
the affine-affine model because with fixed m, n the other ones differ form it by 
(m, n-dependent) c-numbcrs. 

Just as in the classical model, for \n + m\ > \n — m\, i.e., nm > 0, the 
"centrifugal" term 

h 2 (n-m) 2 h 2 (n + m) 2 
cfs "~ 16 Ash 2 1 16Ach 2 f 

is singular repulsive at x — and finite-attractive for |x| — > oo. And then even 
for the purely geodetic incompressible model (V s h = 0) there exist bounded 
states and discrete energy spectrum for E^. For nm < the energy spectrum 
is continuous (scattering states). 
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Finally, let us quote the corresponding formulas for the quantized d'Alembert 
model. Obviously, using the same notation as above we have the following 
expression for the classical kinetic Hamiltonian: 

r d.A = L(p2 + p2) , 1 ™ 2 , 1 " 2 

^ 2I {1+ 2 > + AI(Ql-Q2f + 4/(Ql + Q2)2- 

With fixed values of m, n the problem reduces again to the dynamics of defor- 
mation invariants Q 1 , Q 2 . In the two-polar coordinates the Lebesgue measure 
element is given by 

dl (a; Q\Q 2 ;f3) = P x (Q\ Q 2 ) dQ 1 dQ 2 dad[3, 

where 



Pi = 



(Q 1 ) - (Q 2 ) = \{Q 1 + Q 2 ){Q 1 -Q 2 )\- 



The reduced amplitudes f mn satisfy the eigenequations 

j_jmn jmn rjimn jmn _|_ -y Q 2 ^) f nln J^mn jmn (99) 

with 

fr 2 h 2 m 2 ft 2 n 2 

rpmn rmn y\ £mn j_ rmn , rmn 

' 2/ (/ + 4I(Qi-Q2) 2S + AI{Q^ + Q 2 f ' 

where 

1 d ( d \ 1 d ( d 



The coordinates Q , Q 2 are very badly non-separable even in the very kinetic 
energy expression. There are however other coordinates on the plane of defor- 
mation invariants, much better from this point of view. The simplest ones are 
coordinates Q + , Q~ obtained from Q 1 , Q 2 by the rotation by the angle 7r/4, 

Q+^-^ + Q 2 ), Q-:=-J=(Q 1 -Q 2 ), 

where Q + and Q~ may be expressed in terms of polar and elliptic coordinates 
respectively as 

Q + = r cos ip, Q~=r sintp, 

and 

Q = clip cos A, Q~ = shpsinA. 

In all these variables the Hamilton-Jacobi and Schrddingcr equations with- 
out potential are separable. Obviously, the geodetic d'Alembert model T = 
(I /2)Tr (ip T (p) is completely non-physical. However, the coordinate systems 
(<3 + , Q~), {r,<p), (p, A) enable one to find a class of potentials which are physi- 
cally realistic and at the same time both the Hamilton-Jacobi and Schrddinger 
equations are separable for the corresponding Hamiltonians. 
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The reduced Schrodinger eigenproblem (fM|) with doubly isotropic potentials 
is separable if 

V (Q\Q 2 ) = V+(Q + ) + V-{Q~). 
Namely, we have then H mn = H™ n + H mn , where 

™„ Jl_ ( Q2 , J 3 \ fmn , f ^( m -") 2 , T/ A ,rrm 

~ 27 ^9Q+2 + Q+ 9Q+ J ; + + I 8IQ+ 2 + + > /+ ' 



rrmn = _^ f d 2 J £_\ fmn , ( & ( m + n f , v I fmn 

21 \dQ- 2 Q- dQ- ) \ 8IQ- 2 ' ' 



where 
and 



f mn {Q\Q 2 ) = f+ n {Q + )f mn (Q~), E mn = E™ n + E mn . 



1 + 

Obviously the volume element is given by 

dl (a;Q + ,Q-;(3) =2\Q+\ \Q~\dQ + dQ~ dad(3. 

The doubly isotropic models separable in coordinates (r, tp) are based on poten- 
tials of the form 

V(r,<p) = V r (r) + ^V (p (<p). 
The wave functions are factorized as follows: 

* = e *™<V n/3 / mn (r, 92) = e ma e mf3 R mn (r)$" m (ip) 
and then of course 

= = Hm^, V* = -t* = /in*. 
The reduced Hamiltonian has the form: 

Hmn xrmn , ^ rrmn 
— tl r + -rt , 

r * 

where 

Timn , 3_9\ , „ 

r " 2/ V9r 2 + rarj + 

is as a matter of fact independent of to, n; unlike this, H™" depends explicitly 
on to, n: 

™ _ ti 2 ( d 2 8 \ h 2 m 2 + 2m,ncos(2ip) + n 2 
v -~2l\W 2+ 2Ctg{2lp) d^) + 2l sTn^) + V 

The functions $ mrl ; R mn satisfy eigenequations 

H™(j) mTi — J^mn^mn^ (100) 

H™ + i^™" J iT™ = (101) 
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The ^-equation (|100p is to be solved as first. Then the resulting quantized values 
of E™ n , labelled by an additional quantum number k are to be substituted 
to (jlOip and because of this the labels m, n appear in R although the radial 
operator H™ n , in spite of the used notation, is independent of m,n. There are 
also two additional quantum numbers in R, namely k itself appearing through 
jrjmnk an( j p r0 p er radial quantum number /z, thus, E obtained from (|10ip 
will be denoted by E mnk ^ . 

Let us quote some very interesting model qualitatively compatible with stan- 
dard demands of the macroscopic nonlinear elasticity, 

2k k 2 f 1 Dl + Dl 

V = -a 77r^ + ^r 



r 2 cos(2tp) 2 \D ± D 

In the natural state of elastic equilibrium r = 0, ip = 0. We do not quote more 
complicated and rather non-useful onc-dimcnsional equations for the elliptic 
coordinates p, A. Let us only mention the general shape of separable doubly- 
isotropic potentials 



2 (ch 2 p - cos 2 A) 2 (ch 2 p - cos 2 A) ' 

Finally, we quote a three-parameter family of doubly-isotropic potentials for 
which both the classical and quantum problems are simultaneously separable in 
all the aforementioned coordinate systems: 



-2\ 



2A 2B 



(Q' + Q 2 ) 2 {Q'-Q 2 



+ c((Q 1 f + (Q 2 y 



Here A, B, C are arbitrary constants. It is well-known that the simultaneous 
separability of Hamilton- Jacobi and Schrodinger equations in a few coordinate 
systems has to do with degeneracy and hidden symmetries. 

Two-dimensional models are interesting not only from the philosophical 
point of view of the "Flatland" geometry. They may be practically useful in 
the theory of surfaces of structured bodies and in the dynamics of elongated 
molecules or other structure elements. 
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